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THE  OBLIQUE  WING 
AS  A LIFTING-LINE  PROBLEM  IN 
TRANSONIC  FLOW 


H.K.  Cheng  and  S.Y.  Meng 


SUMMARY 


A transonic- flow  theory  of  thin  oblique  wing  of  high  aspect 
ratio  is  presented,  which  permits  a delineation  of  the  influence  of 
wing  sweep,  center-line  curvature,  and  other  three-dimensional 
(3-D)  effects  on  the  nonlinear  mixed  flow  in  the  framework  of  an 
asymptotic  theory.  The  component  flow  near  the  wing  section  is 
basically  plane  (two-dimensional)  but  nonlinear  and  mixed,  being 
governed  by  equations  consistent  with  the  t ranson i c sma) I -d i s turbance 
approx imat ion . The  work  analyzes  3"D  corrections  to  this  nonlinear 
problem  and  matching  its  solutions  to  that  of  an  outer  flow.  In  the 
(parameter)  domain  of  interest,  the  outer  solutions  correspond  to 
a high  subsonic,  or  a linear  sonic,  outer  flow,  representable  by  a 
Prandt I -Glauert  solution  involving  a swept  (or  curved)  lifting  line 
in  the  leading  approx imat ion . 

Among  the  3*D  effects  is  one  arising  from  the  compressibility 
correction  to  the  velocity  divergence,  absent  in  classical  works. 

More  important  are  the  upwash  corrections  which  include  the  influence 
of  both  the  near  and  far  wakes,  as  well  as  the  local  curvature  of 
the  center  line.  For  straight  oblique  wings,  local  similarities 
exist  in  the  3"D  flow  structure,  permitting  the  reduced  equations 
to  be  solved  once  for  a I I stat  ions . An  analogy  also  exists  between 
the  oblique-wing  problem  and  that  of  a 2-D  transonic  flow  which  is 
weakly  time-dependent;  this  provides  an  alternative  method  of  solv- 
ing numerically  the  inner  airfoil  problem.  Of  theoretical  interest 
is  the  uncovering  of  a compressibility  correction  to  the  outer  flow 
in  the  form  of  a line  source  and  a nonuniformity  of  the  perturbation 
solutions  at  a shock  root  (and  its  treatment). 


A procedure  basad  on  a line  relaxation  method  for  solving  numer- 
ically the  reduced  inner  problem  is  described;  solutions  with  high 
subcritical,  as  well  as  slightly  supercritical,  component  flows  are 
demonstrated.  Comparison  with  corresponding  numerical  solutions 
based  on  f u I 1 -poten t i a I equations  for  oblique  elliptic  wing  shows 
encouraging  agreement.  -r 
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1.  INTRODUCTION 


Use  of  winq  sweep  to  control  the  compressibility  effect  has  lonq  been 
a practice  in  aircraft  design  (see  Busemann  1935,  Jones  1946,  Jones  6 
Cohen  1957,  Kiichermann  1965).  Consider  a high-aspect-ratio  wing  in  a 
supersonic  or  transonic  flow;  there  is  generally  a ranqe  of  wing-sweep 
angle  (for  each  wing  section),  in  which  the  velocity  component  normal  to  the 
center  line  is  sufficiently  low  to  permit  a shock-free  (or  nearly  shock-free) 
component  flow,  yet  is  close  enough  to  the  sonic  speed  for  the  maintenance 
of  a high  lift  coefficient.  In  this  sweep  range  lies  the  potentiality  for 
applying  the  well  documented  supercritical  airfoil  data  (see  Bauer  et  al. 
19/4,  Boerstoel  197**,  Kacpr/ynski  et  al.  1971.  and  Whitcomb  197*0 
to  aircraft  design.  The  flow  field,  in  this  case,  is  necessarily  nonlinear 
and  admits  features  charac ter i z 1 ng  the  mixed  (elliptic-hyperbolic)  flow. 

In  the  present  work,  we  study  the  three-dimensional  (3~D)  structure  of  this 
nonlinear  mixed  flow  by  solving  a perturbation  problem  and  determining  its 
solution  through  matching  with  a lifting-line  solution  to  an  outer  problem. 

The  approach  taken  is  that  of  an  asymptotic  theory  for  high-aspect- 
ratio  planar  wings,  and,  therefore,  shares  a common  base  with  Prandtl's 
(1918)  lifting-line  theory,  as  well  as  the  more  systematic  development  of 
Van  Dyke  (1964)  and  Ashley  f.  Landahl  (1965).  Apart  from  the  nonlinear 
mixed  flow  and  other  compressibility  effects,  the  following  development 
takes  into  account  the  influence  of  the  wing  sweep,  as  well  as  the  cur- 
vature of  the  planform's  center  I i ne( absent  in  the  classical  work.  The 
analysis  was  mot i vated , i n part, by  Jones'  (1971,  1977)  works  on  the 
straight  oblique  wing^which  offers  an  ideal  opportunity  for  extending  the 
lifting-line  concept  and  naturally  becomes  the  focus  of  the  present 
i nves  t i ga  t ion . 

Similar  extensions  have  been  made  earlier  by  the  first  author  in  the 
context  of  unsteady  incompressible  flow  for  a curved  center-line  (Cheng 
1976)  and  of  steady  Incompressible  flow  for  an  oblique  wing  (Cheng  '978a 
and  b) . In  the  latter  work,  agreement  of  the  lifting-line  results  with 

Cheng's  work  involving  a curved  center  line  in  the  steady  incompressible 
case  has  not,  as  yet,  been  published. 


those  from  the  numerical,  wing-panel  and  other  solutions  obtained  for 
high  aspect-ratio  wings  are  found  to  be  consistently  good.  The  present 
study  is  carried  out  in  a framework  equivalent  to  the  transonic  small- 
disturbance  approximation  for  a potential  flow  with  an  aim  of  bringing 
about  a greater  simplicity  in  the  problem  analysis,  which  will  greatly 
reduce  computation  works  and  shed  light  on  the  3"D  structure  of  the  tran- 
sonic mixed  flow.  For  straight  oblique  wings,  the  theory  yields  an  analogy 
with  the  unsteady  flow  of  a two-dimensional  thin  airfoil  in  the  nonlinear 
transonic  regime;  the  formulation  also  admits  a local  similarity  in  the 
3-D  structure,  providing  a further  reduction  of  the  computation  work  for 
a rather  useful  class  of  oblique-wing  geometry. 

An  outline  of  the  analysis  concerning  the  similarity  solutions  of  the 
oblique  wings  has  been  sketched  in  a greatly  condensed  note  (Cheng  6 Meng 
1979)-  This  report  presents  the  theory  i n a ful  1 er  context ; apart  from  the 
crucial  details  in  the  theoretical  development,  the  work  presents  several 
comparisons  of  the  analysis  with  computer  solutions  based  on  ful 1 -potent ial 
equations  for  oblique  wings  at  high-subcr i t ical , as  well  as  slightly- 
supercri t ical , component  flows. + The  theory  developed  is  valid  also  for  a 
curved  centerline  (see  requirements  below)  and  is  expected  to  hold 
in  the  case  of  a conventional  symmetric  swept  wing,  although  without 
further  treatment,  the  solution  is  not  applicable  to  the  vicinity  of  the 
wing  apex. 

For  a straight  unswept  wing,  the  3-D  effect  would  appear  mainly  as 
a local  upwash  correction;  more  specifically,  the  local  flow  field  next 
to  the  wing  section  can  be  taken  as  being  planar  as  in  Prandtl's  original 
theory.  Cook  and  Cole  (1978)  have  given  a rather  complete  elucidation  of 
this  problem;  Cook  (1978)  has  also  studied  the  uniqueness  of  the  solution 
to  the  reduced  problem.  With  the  wing  sweep,  two  important  features  emerge. 
One  is  a pronounced  upwash  induced  by  the  near-wake  vorticities,  which  is 

+ The  report  is  an  expanded  version  of  an  earlier  draft  entitled,  and 
referred  to  in  Cheng  6 Meng  (1979)  as,  "Transonic  Lifting-Line  Theory 
of  Oblique  Wings  and  Other  High-Aspect-Ratio  Planforms". 
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responsible  for  the  pronounced  rolling  moment  charac ter  I s 1 1 c of  an  un- 
twisted oblique  winq  (see  Jones  1977,  also  Cheng  1978a);  the  other  Is  an 
additional  compressibility  correction  arising  from  the  spanwi se  density 
variation  ( ^2.1).  On  account  of  the  latter  feature,  the  flow  field 
next  to  each  wing  section  can  no  longer  be  treated  as  being  planar. 

At  a subsonic  Mach  number,  there  is  a local  spanwi se  component  of  the 
near-wake  vor t Ic I ty,  wh ich  Is  not  negligibly  small  and  must  be  accounted 
for  in  formulating  the  reduced  Inner  problem.  This  feature  does  not  appear 
however  in  the  present  analysis  owing  to  the  relatively  small  sweep  anqle 
required  by  the  nonlinear  domain  (lest  the  inner  problem  becomes  linear).4 

An  interesting  development  in  the  present  asymptotic  theory  is  the 
uncovering  of  a line  source  In  the  higher-order  lifting-line  theory,  un- 
suspected in  previous  work.  The  effect  in  question  arises  from  the 
spanwise  compressibility  correction  mentioned  and  depends  on  the  sectional 
lift  (non  I i near  I y),  as  well  as  the  spanwise  and  chordwlse  distributions  of 
the  wing  thickness  (Appendix  III).  Its  appearance  reminisces  the  nonlinear- 
lift  influence  on  the  transonic  equivalence  rule  (Cheng  £ Hafez  1979)  but 
plays  a lesser  role.  Nevertheless,  this  source  term  is  found  to  be  useful 
in  providing  a description  at  the  far  boundary, needed  in  subsequent 
computations  ( §5)-  In  passing,  we  observe  that  similar  source  effects 
should  also  arise  in  a linear  problem,  although  the  lift  distributions  of 
interest  will  not  be  affected  (to  the  order  of  A,' ) in  such  a case. 

The  theory  admits  shock  waves  and  takes  into  account  the  perturbations 
in  their  strengths  and  locations.  As  well  known  from  Oswatitsch  f.  Zlerep's 
(I960)  work,  there  is  a weak  singularity  at  the  juncture  of  the  shock  and 
the  curved  portion  of  an  impermeable  surface,  where  the  flow  will  undergo 
a reexpansion.  In  the  perturbation  part  of  the  solution,  this  singularity 
manifests  to  a logarithmic  singularity  In  the  pressure  distribution  behind 
the  shock  root.  Althoughthe  validity  of  the  perturbed  solution  In  regions 
removed  from  the  shock  root  could  be  defended, the  breakdown  in  question 
raises  doubts  on  the  ability  of  such  a solution  procedure  for  computing 

+ This  fact,  and  the  two  features  mentioned  above,  have  been  noted  earlier 
by  Cheng  6 Hafez  (1973). 


correctly  the  shock-excursion  distance  (displacement).  Essentially  the 
same  problem  (of  nonuniformity)  undermines  the  shock  analysis  in  the  con- 
text of  an  unsteady  transonic  plane  flow  ( $ . 1 ) where  the  difficulty  is 
compounded  by  other  complexities  in  computations  and  has  not  been  fully 
recognized.  This  problem  with  the  reexpansion  singularity  is  treated  in 

S 2.6. 

In  the  following,  2b  denotes  the  wing  span  (measured  perpendicularly 
to  the  main  flow  direction)  and  the  root  chord  (measured  perpendicularly 
to  the  centerline).  Refer  to  Figure  1.  The  aspect  ratio  is  defined  as 

ib/c,.  We  use  ot  to  characterize  the  wing  camber  and  the  angle  of 
attack  (measured  from  the  zero-lift  angle),  and  will  assume  that  the  wing- 
thickness  ratio  is  also  comparable  to  0 i . The  local  swept  angle  is  A ; 
the  free-stream  Mach  number  is  M..  and  the  Mach  number  of  the  free-stream 
component  flow  is  A'f.cosA  . Three  parameters  will  control  the 

reduced  problems 

Kf  (I  - M„‘  )/*'  , 05  A/cf,  € = '/&*,  ( 1 ) 

where  £ may  be  regarded  as  the  reciprocal  of  a reduced  aspect  ratio,  and 
( is  a rescaled  (typical)  sweep  angle.  The  parameter  is  essentially 
the  transonic  similarity  parameter  written  for  the  component  flow,  with 
the  factor  in  the  denominator  taken  as  the  first  power  of  M„  , as  in 
Murman  S Krupp  (1971);  the  theoretical  reason  supporting  this  preference 
is  given  in  § 2.  The  flow  field  next  to  the  wing  section  and  that  removed 
from  it  are  analyzed  as  two  asymptotically  distinct  (inner  and  outer) 
regions  for  €-+0  with  fixed  K„,  © , and  OL  '*/€  . The  last  quotient  is 
fixed  in  order  to  simplify  the  remainder  estimate.  As  will  become  apparent 
later  in  § 3.  the  relative  magnitude  of  (*)  and  K„  will  also  determine 
whether  the  outer  flow  is  (high)  subsonic,  (near)  sonic,  or  (low)  supersonic. 
Only  the  domain  in  which  K,  > 0 and  K.l  «.-©  0 • corresponding  to  a 

high  subsonic,  or  a linear  sonic,  outer  flow  is  considered. + Since 


Other  regimes  with  Kn  < O , as  well  K0<O  implies  a much  stronger  shock 
extending  far  out  into  the  outer  flow,  and  perhaps  are  less  attractive 
for  aerodynamic  applications. 


the  theory  permits  a curved  center  line,  a smooth  spanwlse  variation 
in  and  (3)  are  also  allowed.  The  result  for  an  unyawed  straight  wing 
in  a high  subsonic  flow  is  recovered  in  the  special  limit  O+o.* 

In  the  next  section,  we  formulate  the  perturbation  problem  for  the 
inner  region  next  to  the  wing  section  and  analyze  the  far-field  behavior, 
the  reexpansion  s i ngu  I ar  i t i es,  and  other  solution  properties.  The  solution 
in  the  outer  region  and  its  matching  with  the  inner  solution  are  presented 
in  Section  3.  The  important  properties  of  the  straight  oblique  wings 
(the  unsteady  analogy,  and  the  local  similarity  in  the  3"0  flow  structure, 
in  particular)  are  brought  out  in  Section  I*.  The  essential  elements  in  our 
computat iona I procedure  are  discussed  in  Section  5.  where  solutions  for 
oblique  wings  with  h i gh-subcr I t i ca I and  s I ight 1 y-supercr i t ica 1 component 
flows  are  demonstrated;  encouraging  comparisons  with  results  from  full 
potential  3"D  solutions  and  from  solutions  via  the  unsteady  analogy  are 
also  shown.  Fuller  descriptions  of  the  computer  work,  with  more  extensive 
results  involving  higher  supercritical  component  flows  and  shocks,  are 
presented  in  separate  works. 


+ A recent  study  by  P.  Cook  presented  at  the  8th  U.S.  National  Congress 

of  Applied  Mechanics,  held  in  Los  Angeles,  June  26-30,  1978,  also  considers 
sweep  effect  in  transonic  flow  corresponding  to  O . Cook's 

paper  is  not  yet  available  for  comparison. 


2.  THE  INNER  AIRFOIL  PROBLEM 


2.1  Governing  equations 

We  consider  a steady  irrotational  flow  of  a calorically  perfect  gas 
with  a uniform  freestream.  The  partial  differential  equation  ( PDE ) for 
the  perturbation  velocity  potential  Ip  can  be  written  in  Cartesian 
coordinates  (x,y ,z)  as 

where  the  remainder  (...)  consists  of  triple  and  quadruple  products  of 
t . e ■ 9 • U.'*'  Pint  ) etc.,  1 is  the  specific-heat  ratio  and  subscripts 
x,  y and  z signify  partial  derivatives. 

For  the  analysis  of  the  flow  field  in  the  neighborhood  of  the  wing 
section,  referred  to  as  the  inner  region  hereafter,  a right-handed  curvi- 
linear orthogonal  coordinate  system  (x1,  y 1 , z1)  will  be  used  (ref.  Fig.  2) 
in  which  x1  and  z’  are  distances  from  the  center  line  (a  reference  curve) 
of  the  wing  planform,  with  z1  measured  along  a normal  from  the  wing  plane. 
The  y 1 -coord inate  is  then  a distance  measured  along  the  (curved)  center 
line.  Writing  A as  the  local  sweep  angle  of  the  center  line,  one  ob- 
serves the  following  identities 


dxt*dyi+dz*=  dK,%f  , ht  - I - x 


(2.2) 


(2J) 

The  PDE  (2.1)  may  then  be  expressed  in  the  curvilinear  coordinates  as 


(I-  m;  )&v.+  <pzt.  - irftaA  h\'<p+r-  ( l-Mlu«A)hmt'Ay.  <P*-  t 

♦ ht  1 £,)  £(  i ♦ ) 4>**  <p*+  U»A  h~t'pM  <py  . + 

+ ( I*  tax  A)  fy'J  + (2. 
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where  U„=U»w<A  and  M#cosA  may  vary  with  y'  and  the  remainder 

....  is  of  the  same  order  as  in  (2.1). 

In  anticipation  of  a local  component  flow  which  will  be  nearly  two- 
dimensional  and  transonic  in  the  x'-z1  plane,  we  employ  the  dimensionless 
var iables 


A 

X ■ 


2*' 


A 

X 


P 


2 j> 


(2.5) 


The  variables  x,  z and  are  identical  to  those  in  the  earlier  works  on 
the  sma I I -di sturbance  theory  of  plane  transonic  flow  (Sprieter  1953;  also 
see  Ashley  and  Landahl  1965),  except  for  the  trivial  differences  of  using 
a half  chord  ^ /2  instead  of  Ct ) and  omitting  fac  tors  with  (it  I ) . Using 
these  variables,  and  eliminating  ( I - M*  ) , tanA  , and  Q /b  through  Kn  , 
(£)  and  € , (2. A)  can  be  reduced  to 


[k,  Pi  - -y  J ♦ $££  = e * ©']  * 


(2.6) 


where  (*)'  is  ; the  terms  omitted  (....)  belong  to  orders  C*  and 

under  the  stipulation  of  a finite  Km  and  (£)  , or  simply  to  the  order 
for  finite  K*  , © and  * %.  The  pressure  perturbation  can  be  computed 

A 

from  as 


(2.7) 


4 

subject  to  an  error  comparable  to  £ . The  retention  of  the  factor  M„ 
in  the  definition  of  Kn  > in  scaling  factors  for  ^ and  * , as  well  as  the 
factor  cosA  in  l/„=  (£ cosA  of  (2.5)  and  (2.7),  are  not  fundamental  from 
the  viewpoint  of  an  asymptotic  theory.  Our  preference  for  these  choices 
will  be  explained  later  (Section  2.3) .+ 

For  the  present  study,  we  consider  impermeable  wing  surfaces  prescribed 
over  the  planform  a(y)  < i < S (?)  as 

*'=  f[«  £*(*'?)  <■  (2.0) 


These  factors  and  cosA  were  not  used  in  the  definition  of  K*  , 
and  % in  our  earlier  note  (Cheng  6 Meng  1979). 
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where  4 and  ft  locate  the  leading  and  trailing  edge,  respectively,  the 
superscripts  ♦ and  - refer  to  the  upper  and  lower  surfaces,  respectively, 

AAA  + A A 

and  Z • Z|  and  I are  all  of  unit  order.  The  functions  1(9)  and  Z*  (?) 
signify,  respectively,  a twist  and  an  upward  deflection.  The  latter  is 

A A 

referred  to  as  "wing  bend"  by  Jones  (1971,  1977).  Both  and  J have 
been  scaled  in  a way  to  allow  the  control  of  the  3'D  effects  (see  below). 
The  boundary  condition  on  the  wing  can  be  transferred  to  the  wing  plane 
(/’  - 0)  as 


(ff)=  £z%«9z ;.«f 


(2.9) 


where  Zj  i s d Z|  / 4y  , and  the  subscript  **refers  to  points  on  the  wing 
plane  bounded  by  0 4 ^ b . The  remainder,  (•••)  including  errors  from  the 

boundary  condition  transfer,  belong  to  orders  or  C*,  same  as  in 

(2.6).  The  pressure  and  the  velocity  component  normal  to  the  trail ing- 
vortex  sheet  is  required  to  be  continuous  across  the  sheet.  These  conditions 
are  again  transferred  to  /’  » 0 to  read,  subject  to  error  comparable  to  £*  , 


ff  $t]  = H$e]I  - 0 


(2.10) 


where  H ]J  signifies  the  difference  across  the  surface  in  question,  and 
the  subscript  TV  refers  to  the  portion  of  the  wing  plane  at  z * 0 
corresponding  to  the  t ra i 1 i ng- Vortex  (Tv)  sheet.  The  first  of  (2.10) 
confirms  that  the  cross-stream  vorticity  of  the  near  wake  is  unimportant 
The  streamwise  vorticity  component  of  the  near  wake  is  nevertheless  very 
significant  (see  (3.12)  later). 

2 . 2 Cha rac  ier  i st_i  c_  and  shock 

The  reduced  PDE  comprised  of  terms  shown  in  (2.6)  has  the  character- 
i s t ic  * z *c(y ■ £ ) satisfying  (C our ant  and  Hilbert,  1965#  pp.  551*556) 

-f  M ^ , A A.  \ 

We  assume  that,  in  approaching  the  lead  mg  edge,  4 /v  t const . lx-a) 
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Kn  - (4*1)  $i  +Z*®  57= -(fj-,)  ; (2.11) 

its  existence  as  a real  surface  requires 

K„  - C1*0  ft  + JJ  < O.  (2.12) 

The  limit  € -*  0 (with  finite  (*)  an d K„  ) determines  a critical  value  of 
the  component- flow  velocity  with 

$a  ~ Kn/(i?0  = U#  (2.13) 


This  identifies  the  transition  boundary  of  the  hyperbolic  and  the  elliptic 
regions,  and  is  also  the  locus  where  the  component  velocity  reaches  the 
sonic  speed  --  a fact  familiar  from  the  transonic  smal 1 -d i sturbance  theory 
for  the  p I ane  flow.  Condition  (2.13),  as  well  as  its  modification  with 
the  0(6  ) correction  in  (2.12),  differs  from  the  critical  condition  for 
the  fuller  3'0  POE  (2.1)  or  its  equivalence  (2.4).  The  transition  boundary 
for  the  fuller  equation  is  simply  the  sonic  condition  for  the  resul tant 
veloc  i ty  in  three  dimensions, + which  derives  from  the  property  of  char- 
acteristic conoids  generated  from  a point.  Whereas  the  cha rac ter i s t i c 
(2.11)  represents  a hyperbolic-wave  surface  generated  by  a line  of 
disturbances  along  the  wing  span,  approaching  a cylindrical  symmetry  in 
the  limit  6-* 0 . 

PDE  (2.6)  and  the  equations  expressing  i rrotational i ty  admit  weak 
solutions  with  surfaces  of  di  scont  inuity  in  V $ , say  x s y , z ) , 

identifiable  with  gasdynamic  shock  waves.  The  shock  conditions  admissible 
to  the  weak  solution  and  consistent  with  Rank i ne-Hugono i t relations,  can 
be  wr i tten  as , at 


(2.14a) 


+ This  critical  condition  is  4>*/u+  , not  reducible  to  (2.13), 

even  in  the  limit  £ -+0  with  /(„  and  © fixed.  Also  note  that,  if  the 
centerline  were  slightly  shifted  to  become  parallel  to  the  characteristic 
surface,  so  that  K„'  s Km  + 26(£)9xe/4/  , (2.11)  and  (2.12)  would  be 
geometrically  invarient. 
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(2. 19b) 


n $ JJ  * 0,  (2.1  9b) 

where,  as  before,  < > and  a b stands  for  the  arithmetical  mean  and  the 

jump  of  the  quantity  across  the  surface  in  question,  respectively.  The 
* 

continuity  of  ^ ( 2 . 1 4b ) , implies  the  continuity  of  the  tanqential 
velocities, 


Ifr I •'  ff  & A-  Hfrfl  - -/ : §7  57 


(2.15) 


2 , 3 On  choos i n g t h e form  for  Kw  and  retaining  cos  A i n I7„  . 

The  use  of  H for  the  2-D  transonic  small- 

disturbance  solution  has  been  proposed  by  Murman  and  Cole  (1971),  also 
Murman  and  Krupp  (1971),  which  yields  consistently  good  agreement  with 
exact  solutions  for  transonic  flows  past  shock-free  supercritical  airfoils 
(Korn  1969).  Apart  f rom  practical  considerations,  the  very  similar  para- 
meter^in  the  present  paper  w i t h ^ rep  1 ac  i ng  i s adopted, 

however,  for  additional  theoretical  reasons. 

A 

We  observe  that  terms  deleted  from  the  coefficient  of  <p  $ £ in 
(2.6)  will  give  an  incorrect  description  of  characteristic  slopes  near  the 
transition  boundary  and  may  lead  to  a nonuniformity  of  the  inner  solution 
(at  some  stage  of  the  higher-order  expansion).  This  deficiency  can  be 
improved  through  a correction  to  the  critical  speed  for  the  component  flow, 
hence,  to  the  K„  in  (2  .6).  The  latter  can  be  explicitly  determined 
independently  of  the  solution  (Appendix  l).  The  result  takes  the  form 
of  a slightly  modified  definition  for  the  similarity  parameter,  namely 

Kn'  * (t-nj  )/<*V*M*  , Ulzftit, )/}(«/)  (2.16) 

+ See  Bailey  6 Ballhaus  (1973)  for  similar  consideration  in  3*0  cases. 

Among  those  deleted  from  the  coefficient  the  major  terms  are 


'l«-0Kn<Pi  -lf(zi-i)  tf]  . 


The  value  of  u>  ranges  from  5/6  to  23/24  for  1 = 1 to  5/3,  and,  therefore 
has  been  taken  as  un i ty  in  our  definition  of  Kn. 

A similar  transonic  parameter,  though  in  a more  complicated  form^has 
been  pointed  out  by  Sirovich  and  Ho  (1978),  who  considered  simple-wave 
solutions  in  a purely  supercritical  case,  for  which  a nonuniformity  at  the 
transition  boundary  can  be  readily  identified.  In  a study  by  Klunker  and 
Newman  (1974),  a term  - is  added  to  the  coefficient  of  j>KA  of  POE 

(2.1)  with  an  aim  to  better  approximate  the  critical  speed,  the  result  may 
be  shown  to  be  equivalent  to  a modification  of  the  similarity  parameter  in 
the  manner  of  (2.16)  for  a straight  unyawed  wing. 

The  component  of  the  free-stream  velocity  U„=  T£cos A which  appears 
in  (2.5)  and  (2.7)  defining^  and  ty  could  be  replaced  simply  by  caus- 
ing errors  no  larger  than  0(  € ^)  or  0(  oi ^ ),  since  A — = 0 ( ot  ) 

has  been  assumed.  However,  its  retention  in  the  i nner  solution  has  a 
definite  advantage  for  enlarging  (slightly)  the  applicability  range  of  the 
sweep  angle.  This  follows  from  an  examination  of  terms  associated  with  A 
in  the  remainder  of  (2.6)  - (2.10),  comparing  them  to  the  3*D  effects  of 
interest.  The  latter  belong  to  0(6  0 ) . The  errors  resulting  from  re- 
placi ng  Vn  wi th  U.  in  (2.5)  and  (2.7)  are  proportional  to  _A  • The 
remainders  related  to  A in  PDE  (2.6)  and  in  (2.7)>as  well  as  (2.10),  fall 
into  two  types,  one  is  0(  e©A*)  and  the  other  0(  A/*,).  The  relative 
magnitude  of  the  four  groups  in  question  are ^there fore , i n the  proportion: 

£0  : A2  = €0A*  : /)/*;'  = I : ■■  e1®’  ■ € , t2-1?) 

i/| 

where  the  assumption  of  € = 0(  eL  ) is  made  as  before  and  the  order  of 
magnitude  symbols  have  been  omitted.  From  (2.17),  it  is  clear  that  if 
u.  were  replaced  by  Vgg  , errors  in  the  solution  would  be  of  an  order  € 
higher  than  the  3"D  effects  analyzed  only  when  (*^=0(l).  Now,  the 
retention  of  the  cosA  factor  in  the  definition  of  and  Cp  eliminates 
the  need  of  considering  the  second  group.  Therefore,  the  error  of  the 
solution  can  remain  an  order  € higher  than  the  3"0  effects  analyzed  as 
long  as 


0 = 0(cw;=0(»“). 


(2.18) 
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or  mii<j  I lor,  according  to  (.’.I/) 


2 . !)  I’erturhnt  Ion  ot  tln>  lomponenl  t low 
I (|u.it  Ions  (2.b),  (2.‘l)  and  (2.10)  pr*  i m 1 1 tin  expansion  ot  (f  in  • 


= <P,  ♦ * 4>, 


♦ • 


(2.  H) 


will'll'  .1  Wl'.lk  (loi)ill  Ithlllll)  llt'pt'llllf  lie  I'  ot  , I'll  ( Is  anticipated.  III!' 

A 

IcMiliiu)  term  P is  .*  0 solution  to  thr  vomuonent  flow  »it  tho  spun 

9 

A 

st.it  ion  y s,itistylni|  tin*  I'OI 

[k.  - i] & t ' fr-  i = 0 

■iii.I  thi'  conditions  ini  tin'  s .isls 


(i?0  - 1 4 ^T;  °i  (2-*>n 

wlii'ii'  tin'  suhsc  l I p t s " ^ 11  ,in.t  "IV"  .,iii\  sonu'  ini'. in  I iw|x  .is  before.  I In'  io 

/\ 

efficient  is  governed  I'V  thi’  I I no.  n nonlioiiKnii'iieous  I’OI  In  vat  I .it*  les 
x .iii.I  / (ol't  .1  i nt'il  .-it  1 1 • i 1*1)11.1 1 I ni|  1 1*  i ms  p I o)»oi  t I oiid  I to  the  (list  | 'owe  i 
ot  « ) 


(2.22) 

.iml  tlie  i on.l  It  Ions  on  the  x axis  (ol't  .i  I ne.f  In  .i  slmll.ii  mantlet  ) 

; Hit (.1-1,^. I'  0,  <*■»> 

where  y appears  essentially  as  a I'lii.imetet  We  note  .i  I so  th.it  the 


It  Is  rather  well  known  th.it  the  t'OI  .joveinlmi  the  ttanxonlc  component 
flow  around  a thin  yawed  wlnu  ot  Infinite  as)»e( t latlo  cannot  he 
recovered  completely  from  I’OI  (2.1)  without  elthn  let.ilnlni)  or 
invokiiu)  the  approx  I mat  Ion  n'sj\  It  I We,  there  foie,  may  expei  t the 
present  work  to  he  slightly  moi  «•  militate  than  those  based  on  the  standard 
VO  transonic  smal  1 -dl  sturhatu  e notation. 


I ^ 


A Ai  ^ ^ ^ ^ . 

solutions  p(  and  shore  the  same  characteristics  system  XsKg  (/'*) 
given  by 


=-(#/ 


(2 . 2M 


Thus,  the  transition  boundary  separating  the  elliptic  and  hyperbolic 

» a A a 

regions  is  K —(4*)Ar4>  ~ 0 for  both  and  <P  . 

fl  fji  ' • O' 

The  equation  describing  the  shock  is  expanded  as 


^ A A \ / A A V A ^ ^ \ J 

X = X (y,  z)  = x0  (y,*;  t e K,  o,r;  + •••, 


(2.25) 


/S® 


allowing  again  a weak  (logarithmic)  dependence  of  X , on  € . The  Hugoniot 

/A  Ap  /\  y\0 

relations  (2.1*0  can  be  transferred  from  X = X to  X-*0  . yielding, 
after  equating  terms  of  equal  algebraic  powers  of  € 

*.-«•<#.>  = -6£)\ 


2 z 2 Z , 


II  £ * J = 0 


(2.27)' 


As  noted  earlier,  the  subcritical  flow  right  behind  the  shock  root 
over  a curved  (wing)  surface  has  a weak  singularity  in  the  surface  value 
of  Alt  , 1 i ke  (x ) , this  will  then  give  a logarithmic 
s i ngul ar i ty  in  in  view  of  the  shock  condition  (2.27).  This  nonuni- 

formity is  treated  in  §2.6. 


For  flow-field  computation,  i t may  be  more  convenient  to  apply  the 
first  of  (2.26)  and  (2.27)  with  (£"_)£  X$,°)£  and  (if)  i el imi nated 
through  the  second  of  (2.26)  and  (2.27)  » e.g.  K(f.)£  +Wk(P.hH 


r 0 


1 k 


l 


2 . 5 Behavior a t large  d i stance 


Since  9 1#/?*  vanishes  as  | <•  | £ J A*  ♦ K„  £'  — ► co  PDE  (2.20)  approaches 
the  Laplace  equation  in  * and  I the  leading  term  in  an  asymptotic 

j'  A 

expansion  for  ^ at  large  |£|  may  therefore  be  represented  by  a vortex 
potential  In  the  variable  £ a x + t 7 • This  and  the  successive  terms 

A ^ 

in  the  development  of  <P0  for £/»  / can  be  inferred  from  POE  (2 . 20 ) under 
the  requirement  that  and  are  continuous  everywhere  Including 

the  x - axis  ( removed  from  the  wing  section) 

$ ^ ,4  wh) ♦ * %0i  1*1  * <£]£.  * 


+ [$;t  r 3: fc£Jt‘r 


(2.28) 


A /N 

where  II  4o  JJ  yg  ' s the  circulation  determined  by  the  potential  jump  at  the 
trailing  edqe(TE);the  remainder  (...)  is  at  most  of  the  order  (I„£|)  f* 

The  coefficients  Dor  end  {5,1-  may  be  taken  as  the  real  and  imaginary 
parts  of  a complex  doublet  strength  Da—  + t D,  ; they  are  functional 

of  and  (#j)‘  , and  can  be  explicitly  computed  in  the 

forms  of  line  and  area  integral s, unaffected  by  the  presence  of  an  imbedded 
shock  (Appendix  II ) . A form  comparable  to  (2.28)  has  been  given  earlier 
by  Cole  (1975). + for  numerical  analysis,  these  explicit  forms  for  the 
doublet  strengths  are  unnecessary,  since  they  can  be  determined  along  with 
solutions  at  and  near  the  far  boundary,  with  the  help  of  a least-square 
method  ( §5),  as  has  been  in  Hafez  L Cheng  ( 1977a  f.  b)  . 

The  PDE  (2.22)  governing  £ approaches  a Poisson  equation  at  larqe 
1*1  , which  gives  the  two  leading  terms  In  the  expansion  of  . These 
and  the  successive  terms  fulfilling  the  continuity  requirements  on 
and  can  be  determined,  for  l£l  » I . as 


The  term  proportional  toIT^'*  in  Eq.  (2.28)  differs  from  that  given  by 
Cole  (1975)-  The  angle  of  the  arc  tanqent  is  restricted  to  the  range 

i-"/t  ,*/*!. 
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$ » @£  [ (4 . 1 a,  fi')i  !„&  * ee'ffl  41 7 . 

' L'®/  2 KM0/2n  2 2H  ffl*J 

- i;R  }*£;:.  ft, . £ u/fi  ♦ 

/Kjfl'L  \ff  l *,  @ ) °o  i *r  g-  ; 0.  J 


/ A * 


^ * 

*-£;  ^"’(rkt) + is9*£ ] f ?(?)  ♦ •••  (2.29) 

where  'sfr£(Urt  '.  (£*  and  <£' correspond , respectively,  to  the  induced 
upwash  and  a component-Mach-number  correction  to  be  determined  later 
through  matching;  the  parameter  K0mK^  “ 0*  can  be  identified  as  the  tran- 
sonic similarity  parameter  at  zero  sweep  ( / - ) / 0 M—  • The 

remainder  (...)  vanishes  as  |^|  — i*  00^  and  the  symbol  designates  terms 
arising  from  nonlinear  corrections  in  the  far  field  of  orders  ( )* 

U|£|  and  unity 

^ ~ 57  f*  Jj  * ?§)  K [*  (t’Wf  * 1 (-jp  - 1 

or* 


(2.30) 


where  S5X+tJicmz  and  £ zl  x -l]Km -z  The  appearance  of  the  source  term 

in  the  far  field  with  a strength ^ £0  should  not  be  too  surprising,  since 


Among  the  remainders  are  terms  of  order  ( U l£l  )*  and  ( i*  |£)  )^/|£  | 
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(2.22),  cast  into  a Poisson  equation 


4 $ = f10»7 ' 


(2.31) 


indicates  clearly  a source  distribution  which  gives  a nonvanishing 
total  vo I ume  flux  to  the  component  flow. 

The  influence  of  this  term  on  the  outer  flow  is  quite  similar  to  the 
one  in  the  transonic  equivalence  rule  involving  lift  (Cheng  6 Hafez,  1973) 
although  the  line-source  effect  has  a considerably  greater  influence  therein. 

A 

The  source  strength  Q ^ is  completely  determined  by  the  far-field  property 
of  4>0  . more  specifically.  D'  and  0f.)!  . In  cases  where  ( IFof  can 
be  neglected,  the  area  integral  of  the  RHS  (2.31)  over  K and  rj  can 
be  correctly  equated  to  the  total  volume  flux  and  determined  as 


Q=  — {( 

K*  'I  iQ  K*  7 


‘ K.  O j . (2.32) 

confirming  the  existence  of  a *0  ( uni  ess  p'  or  the  wing  thickness 
vani shes) . For  (2.29)  shows  that  and  , in  this  case 

cannot  be  evaluated  simply  from  the  area  integral  of  RHS  (2.22),  which  may 

be  unbounded;  the  quantity  of  interest  may  nevertheless  be  determined  by 

£ /* 

constructing  first  f unc  t ions  and  4>e  wh  i ch  satisfy  two  Poisson  equations 

and  agree  with  <^#  and  <£(  , respectively,  in  the  far  field  in  accordance 

with  (2.28)  and  (2 . 29);  app I y i ng  next  a result  of  Green's  theorem  to  the 

A ^ A ^ 

difference  of  ( t • <f, ) and  ( ).  Using  (2.11*),  the  result  can 

be  shown  to  be  unaffected  by  the  presence  of  shocks.  The  explicit  relation 
connecting  Q,  to  Q'  and  Of.)2  will  be  given  later  in  £ A in  the  case  of 
oblique  wings  under  local  similarity;  its  derivation  is  described  in 
Append i x III. 

The  demonstration  in  (2.32)  indicates  that  similar  line-source  effect 
will  arise  in  the  corresponding  lifting-line  theory  in  a linear  subsonic 
flow,  although  the  lifting  and  the  thickness  problems  in  this  case  can  be 
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separately  analyzed  on  account  of  the  linearity. 

A A 

It  Is  assumed  that  the  solutions  ^and^are  uniquely  determined  after 
fulfilling  the  Kutta  condition  at  the  trailing  edge,  the  conditions  (2.21), 
(2.24) , (2.27),  and  the  far-field  behavior  represented  by  the  leading 
groups  of  terms  on  RHS  (2.23)  and  RHS  (2.29).  The  latter  is 


- <r;m  t . 


(2.33) 


The  first  of  (2.33)  is  equivalent  to  the  vanishing  of  the  qradient  of  fig 
and  the  second  requires  specifications  of  the  upwash  and  component  Mach- 
number  corrections. 


2.6  Treatment  of  the  jion unjjo rm i t y at  the  s hoc k root 

Recall  that  the  logarithmic  singularity  in(^>#)<  is  a result  of  the 
breakdown  of  the  analytic  continuation  from  the  perturbed  to 

the  unperturbed  shock  boundaries  , cf.  (2.27)  . I t is  apparent  that  terms 

1 i ke  |(*  *)*<  J a and  x * < f } will  arise  in  the  corresponding  shock 
relations  in  the  next  order.  This  and  similar  reasons  lead  to  the  ex- 
pectation that^^-fis  singular  like  1/x'with  f's  x - X®  , ond  (^j)^like 
’/S'*.  Thus,  a nonuniformity  of  the  expansion  (2.19)  is  identified  at 

S-t‘.  = 0(e)  ■ 

Instead  of  carrying  through  the  formalism  of  expansion  in  a new  inner 
variable  for  this  (inner  of  the  inner)  region,  we  treat  the  nonun i formi ty 
directly  in  the  framework  of  the  composite  system  (2.6),  (2.9),  (2.10) 
and  (2.14a  6 b) . We  shall  first  establish  below  that  the  reexpansion 
singularity  in  (th  , written  in  coordinates  fixed  to  the  shock  root  is 
unaffected  by  the  3_0  corrections.  Using  this  as  a valid  flow  description 
near  the  shock  root,  it  becomes  possible  to  determine  the  3“D  corrections 
to  the  shock  jump,  the  shock-excursion  distance,  as  well  as  another  para- 
meter of  the  reexpansion  singularity,  from  matching  the  perturbation  of 
this  solution  with  the  logarithmically  singular  field  of  ft  > # . 
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For  this  purpose,  we  identify  the  location  of  the  shock  root  by 
X * • we  use  the  subscript  "sr"  to  denote  values  at  x — *jr  on 

the  wing,  and  the  superscripts  and  "+"  to  denote  the  pre-shock  and 
post-shock  cond i t ion,  respect ivel y . For  convenience,  we  introduce 

/3,=|x;-c*'><s;|*4 

Af  A a.  . 

where  Utr  is  the  post-shock  value  of  at  the  shock  root.”  In  addition, 

we  I e t 

D = ([©$*©']  , 

Zws  &Z  * *<®Zt  *tl  . 

PDE  (2.6)  admits  a development  near  the  shock  root  (on  the  subcritical, 
downstream  side)  consistent  with  the  impermeable  wing  condition  (2.9), 

$ = L+ 

C = d C Aa 

where  "R.P.”  stands  for  the  real  part,  and  the  remainder  (...)  is  comparable 

to  0(15,1').  The  latter  corresponds  to  a relative  error  comparable  to  0(0 

in  the  inner  region  c ,=0(«).  The  wing  boundary  condition  also  requires 

' 

that 

I p.(A)-0,  I P(B)  = -(?*),, Afi,,  (2-37) 

where  I.P."  refers  to  the  imaginary  part,  and  stands  for  3L  . The 
real  parts  of./)  together  with  the  shock  geometry  near  the  root  e,  = fA?.> 
can  be  determined  consistently  from  the  two  Hugoniot  relations  (2.114a  6 b) 
and  a regular  upstream  supercritical  solution  for  small  £ and  . The 

, - - -- 

*f  Note  that 
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f i na I resu I t for/\  is 


J)b  (Z:\r/f}, 


(2.38) 


The  real  part  of  Q can  be  determined  only  through  matching  with  the 
solution  in  the  region  removed  from  the  shock  root.  The  3*D  effects  can 
enter  (2.36)  only  through  R.P.  (B)  and  through  the  known  constant  C ■ 
Whereas  the  description  (2.36)  differ  from  the  original  (Oswatitsch  £ 
Zierep  I960)  form  in  the  addition  of  a nonanalytic  tern)  proportional  to 

• A 

K.r the  field  of  the  pressure  perturbat ion^or  <P  £ , is  nevertheless 
analytic  and  of  the  same  original  form.  The  result  can  be  expressed, 
wr  i t i ng  if  as  u , as 


A A 

u - Utr 


HuU, 


'Sr 


~*-p£w< , 


(2.39 a) 


? = {»«'?;  J . 39b> 

Ff  I riiy  " ' 

The  parameter  K,  is  an  arbitrary  constant  containing  the  undetermined  R.P. 
(B)  as  well  as  a contribution  from  D u*r  , The  corresponding  form  of  the 


shock  boundary  is 


~ v2[ wlh(£)  * / J , 


? *x 


St 


(2.A0a) 


(2 .AOb) 


The  above  description  confirms  that  the  analytic  form  of  the  pressure 
A 

perturbation  on  pf  near  the  shock  root  is  unaffected  by  the  3~P  effects. 

As  a matter  of  fact,  this  conclusion  holds  even  if  the  3*0  effect  through 

D u/r  was  of  order  unity  (cf.  (2.35)).  This  observation  will  have  relevance 

in  the  analagous  problem  of  unsteady  transonic  problems  (see  $ A . 1 below). 

The  result  (2.39)  indicates  that  any  3"D  influence  on  this  (inner) 

structure  can  result  only  through  changes  in  the  three  parameters:  irtfij  . 

* *r 
^J^fand  K, . Conversely,  the  structure  determined  through  matching  with 


! i 
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the  perturbed  solution  in  region  exterior  to  this  neighborhood  should 
yield  knowledge  on  the  3~D  corrections  to  aft! 


Sr 


xp 

*tr 


and 


Introduce 


A A A 

5u  I u-u0 


i */,«  *tr  - (*elr,  SK,M  K,  -(*, ),  . 


(2.M) 


The  subscript  "0"  refers  to  values  based  on  f . At  sufficiently  large 
Iff.  we  can  expand  (2.39)  for  small  $ U > £ */r  ar|d  SK,,  taking  into 

account  of  their  effects  on  ^ and  r?  . The  final  result  yields  a logarith- 
mically  singular  behavior  for  the  correction  to  u a V£ 


A + cl  + * d^aryg)* 

+ e%  * e,  T}  + ...  , 


(2.1.2) 


with 


$ = - #(z;)„  (Vf'lutef's* 


A P 
Sr 


*'-tl 

d, =  dl  - -fr  6lIUir/tUo1lsr  , ^1=  w[&l*  ^m]  **Sr  , 

+ Ci-  2- *c.)  [ <j,  -^ssuHsr/a^Vtr]  , 

e.  = <*t  - . 

The  result  (2.42)  corresponds  to  the  inner  solution  (2.39)  in  the 
outer  limit,  and  should  match  with  the  logar i thmi ca 1 I y singular  behavior 
of  the  solution  near  * a K a O.  Since  the  coefficients  etc . 

are  completely  determined  once  the  exter i or  po  and  ^ problems  (posed  in 
ff  2 . 4 and  2.5)  are  solved,  the  matching  then  will  determine  the  3"D 
corrections  to  the  shock-root  location  from  the  velocity  jump  at  the 

root  from  J.  , and  the  parameter  K,  from  © . The  matching  of  the  coefficients 
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d(  and  furnishes  additional  checks  on  the  theory's  internal  consistency. 

In  passing,  we  point  out  that  the  singularities  at  the  leading  and 
trailing  edges,  and  the  related  problems  of  nonun i form i t y , have  not  been 
considered.  It  is  well  known  that  for  a ^JTuconst.  (£-£)*,  the  2-P 
surface  pressure  is  singular  like  (x-f  ) according  to  the  transonic  small- 
disturbance  theory  (Nonweiler  llIS8,  Cole  l<171>).  Keyfltz,  Melnik  and 
Gr<»ssnK»n  (19/8)  are  able  to  show  that  as  long  as  oi  does  not  far  exceed 
the  thickness  ratioCfi.e.  <*/t  = 0(1).  the  symmetric  behavior  with  the  - 
exponent  remains  to  dominate  the  surface  pressure  near  such  a leading  edge. 
One  may  then  argue  that  nonuniformity  cannot  arise  from  f-D  corrections, 
because  the  flow  asymmetry  will  not  alter  the  -*/%  exponent,  and  also 
because  the  spanwise  compressibility  correction  in  the  right  of  PDE  (2.22) 
is,  at  worst,  comparable  to  which  can  be  no  more  singular  than 

terms  on  the  left  of  (2.22)  involving 


Preliminary  computat ional  results  based  on  a shock- f i t t i ng  procedure  to 
determine  and  In  a supercritical  case  Indicates  that  such  a 
procedure  for  determining  the  corrections  for  ffuDjr  , and  K,  Is 

feas i bie. 
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3.  OUTER  SOLUTION  AND  MATCHING 

3 . 1_  The  ou ter  solution 

In  the  flow  domain  of  interest  (K„>0 , K = K„*0}O)the  basic  outer 
flow  is  described  by  the  linear,  elliptic  Prandt 1 -Glauert  equation.  For 
subsequent  analysis,  we  shall  employ  a set  of  dimensionless  variables 

xsVBb.  4>m<t>/d?UmC.,  (3.D 

w i t h _ i r . 

5 , oL , = OC/^  ( 

and  assume  an  expansion  for  a dimensionless  perturbation  potential 

(3-2) 

where  a weak  (logarithmic)  dependence  of  4*,  on  £ is  anticipated.  The  term 
{ accounts  for  the  higher  aspect-ratio  effect  as  well  as  the  nonlinear 
correction  to  the  outer  solution.  PDE  (2.1)  yields  equations  governing 
>ind  41, 


(£  * ♦ il  \ A = o 

\ 7il  if1  >*'  / • 

(L  > 11  = JIL  ±{i  -)* 


In  the  far  , i eld,  the  gradients  of  ^ and  are  required  to  vanish, 

and  on  the  wing  trace  behind  the  center  line,  the  continuity  of  the  pressure 
and  normal  velocity  across  the  t ra i I i ng- vortex  sheet  dictates  that  the 
partial  x and  * derivatives  of  and  4>  be  continuous  there. 

The  <pt  is  the  velocity  potential  of  a lifting  line,  representing  the 
bound  and  trailing  vortex  system  associated  with  the  center  line  *=Xc(y)  + 


f'.L.f  £W—.rlt 

4 l,  (r-ty-  ■ * L «, 

where  R,  e [ (i-  ic if, )j**(9  -?■)*♦  * 1 j'? 


S-5c  </. ) 


Jd\ 


+ The  factor  * in  (3.M  was  left  out  in  Eq.  OM  of  Cheng  1978a, in  Eq . (5.0 
of  Cheng  1978b,  as  well  as  In  Eq.  (12), Cheng  £ Meng  1979  ~ an  unforgivable 
transcribing  error  called  to  our  attention  by  Dr.  T.  Evans. 
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In  approaching  the  center  line,  i.e.  ~ 1 * 0 . ^ "*& , the 

integrand  shown  in  (3.3)  becomes  non integrabl e at  /,  = X . By  subtracting 
a suitable  function,  say  g,  from  the  integrand,  the  resultant  integral  may 
yield  a limit  as  ? and  *•  vanish.  This,  together  with  the  quadrature  of 
g,  then  gives  the  behavior  of  in  the  inner  limit  sought.  This  procedure 
has  been  used  in  the  analysis  for  straight  oblique  wings  (Cheng  1978a  6 b) . 
Among  the  indefinite  integrals  in  the  quadrature  mentioned,  the  following 
two  are  bas  i c 


= (C5) 


wht  rt  3nd  ^ ii r c 

50*  - wg/f/  ♦*»*),  (V*  *')/( /♦"*),  ] 

I (3.6a) 


•¥ 

Among  the  typical  definite  integrals  resulting  from  the  quadrature  are 
(omitting  all  the  bars)5 

l.jr  I-/ 

I Z f !+ ($-»*•) /JTim*  R”]  du  = firj  ^ I ^2  [ *•*(*/$*&  t)  f 

-i-/  -'-r 

( /*»*)  4 m ) //-  y l''+(fi-tmk  •**•)/»/ 1 ~'J 

f /7**»  [f  l*  (\ /(*''*))  ~ ( 

-l-f  ^ 

(e  ttfifm1 1)  + **  [t*» '(*/(, -J))  + * 


R a i b.u  t Z, 


\ ■ P ♦ ‘ Q 


(3.6b) 


P"C.  t s,u  - R p /T.  R,  Q * r(B0  + u)  - I P/?'- R 


(3.6c) 


The  equation  xzx((y)  defines  the  centerline  in  the  Prand  1 t -G 1 auer t 
variables. 

To  arrive  at  a suitable  form  for  matching  with  the  inner  solution, 
the  result  for  sma II  £ — x - Xt  ( y ) and  f is  first  trans  formed  to  a 
system  of  curvillinear  orthogonal  coordinates  ( H , y e with 
Z — Z , and  the  curved  /-axis  identified  with  x = XcCjr)  , also 

vr,*  -at  dA  _ ,, 


,_i  J -r 2 , / -'A  A \,-1  dA  fh  (y  1 

d?  * df  -dx  t (l  - X dA_  \ dy  — - . .. 

V dy>  ' ' dY’  (/#*»/* 


This  is  followed  by  a t ransformat ion  to  the  inner  variables  X ,z  and  y' t 


observing  that 


CM  = l&A  = IK  Sin  A = I*£  = §L 
Jl ^ IK*  ’ K 

= L dJL  = A i_f  «*djZi  = K 

B • dy‘  (I -HA)'1  ' dy  0 " ^ 


(3.3a) 


and  that;near  the  inner  limit 


X = — ! — JL  . M.1  eanVl  dA  *1  + 

r'zj  B'P^'y  * ;ir,  + 


(3.8b) 


The  final  result  for  in  the  inner  limit  is 


where  the  remainder  is  comparable  to 


«/<?!*,  ' and 


£r' 

K, 

+ 


(„  I !_L!  I t 2]  + ir.il  * UM|< 

(coiVi  I J <*?'  I •*?  I -Sin At 

+ f } [ I - ty(y,-?)]  jy.  - 

i,  /.  - y 

{ r.(»)[  . - 


- 5/r.i  J5  „(y- -y)J-i^C^)/7,-y|j 


(3-9b) 


noting  that  Jy'/dy  = sec  A ,Ko=IC-02  , and  that  we  have  replaced 
cos  A (notCcS/1  ) by  unity  and  /'by  / , throughout  (3-9)  except  in  the 
leading  term  (cf.  § 2.3).  ^ ^ 

The  outer  expansion  of  the  i nner  solution  ( ^ + € 4>,  + * ” ) , accord i ng 
to  (2.28)  and  (2.29),  share  a domain  of  validity  with  the  inner  expansion 
of  the  outer  solution  ( <£o  ♦ — ) in  f «/<$  | « t'1  where  ma  tch  i ng  is  achieved 
after  identifying 


no) = r. «f))  = cosa  r„  (/ );  c.  tx)  = o,i 
<£  - (*,  j v‘  [vr-(«»f.)"i>*aNil')]: 


(3.io) 


The  induced  upwash  and,  hence,  the  inner  solution  is  completely  determined. 
Terms  in^  + e^not  matched  at  this  stage  are  those  from  the  inner  solutions 
belonging  to  order 


Iff' 


«/i»  ?r , (bn 


e 


(3-11) 


These  should  find  their  counterparts  in  4>,  which  includes  the  nonlinear 
corrections  to  the  outer  solution. 


+ For  straight  oblique  wings,  is  precisely  the  £ in  Cheng  and 

Meng  (1979).  * 


h 

i i 

H 

? i 

s 

M 

■j 

. . 
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The  upwash  function 


v; 


is  dominated  by 


£ 

*y 


r ( SinA  I')  H£l , 


(3.12) 


the  part  proportioned  to  dro/j?  represents  the  important  influence  of  the 
near  wake,  the  part  proportioned  to  represents  the  downwash  induced  by 
the  bound  vorticities.  Those  together  with  other  competing  contributions, 
lead  generally  to  a maximum  induced  upwash  near  the  wing  tip  on  the  aft 
wing  panel  of  an  obligue  wing  and  a maximum  induced  downwash  near  the  t i p 
on  a swept- forward  panel,  as  noted  earlier  by  Cheng  (1970b,  p.!7).+ 

In  the  very  special  limit  corresponding  to  an  upyawed  straight  wing 
(©♦O).  the  upwash  correction  determined  from  (3.9)  and  (3.10)  agrees 
with  the  linear  subsonic  results  of  Jones  f,  Cohen  (1997).  Interestingly, 
for  a symmetric  swept  wing,  i.e.  = mly  I , the  last  integral  of 

\J°°  diverges  in  approaching  the  wing  apex  like 

1 ' 4n  y ■ 

This  agrees  with  the  downwash  induced  by  the  bound  vortices  on  the  opposite 
wing  panel  computed  according  to  the  Biot-Savart  law.  It  exhibits  clearly 
the  extent  of  nonuniformity  at  the  apex 

y = 0(©  K«V,0.  (3.13h) 


The  remainders  in  the  matching  listed  in  (3.11)  contains  contributions 
from  the  source  terms  as  well  as  nonlinear  corrections,  which  are  of  theo- 
retical interest.  To  ascertain  more  clearly  their  influence  on  the  outer 
flow,  we  shall  investigate  in  the  following  the  correction  to  the  lifting- 
line solution  4 \ • and  examine  its  matching  with  the  inner  solution.  To 
avoid  tedious  details,  only  t lie  case  with  a straight  center  line  will  be 
cons i dered . 


Same 
w i ng 

dA/dy ' 


can  be 
With 


said  of  the  upwash  field  near  the  tip  of  a symmetric  swept 
a suitably  chosen  distribution  of  the  center- line  curvature 
, such  an  upwash  effect  can  be  mitigated  (Cheng  197b). 


27 


3 • 3 Higher-order  outer  solution  and  matching 

The  second  PDE  of  (3-3)  governing  is  a 3-D  Poisson  equation. 
In  the  curvilinear  coord  I na tes ( x' , y ' and  z'  ) of  (3-7).  this  equation 
becomes  ( with  ^ i=  • - i'dX/d y ' ) 


( hi  4 b‘ ) i * 1 «*«($) . 


(3.1^) 


We  shall  confine  the  analysis  to  the  caseQ's  0 ( t he  swept  wings  with  a 
straight  center  line)  for  which  the  term  with  dti/*}  on  the  IMS  (3.1M  is 
absent . Writ i ng 


(3.9)  ijives  for  sma  I 


*5  = x'  + i £’(  7'-i'2' 


( 3 . 1 ‘i ) 


rJ-  ~ -£)«■•• 

witii  the  remainder  (•••)  vanishing  as  i<r 'b0*.  Thus  ( 3 . 1 /* ) may  take  (tie 
form  for  small  in  "• 

w,  * ~ f - 

- c*5\l  rr'J  ro  d P,/4y'  (~p  * ^ )J+  (3>,6) 

where  the  remainder  (•••)  includes  contribution  from  / 3 / ' J 

Integrating  with  respect  to  £ and  , (3- lb)  yields  the  asymptotic  be- 
havior of  the  particular  solution  for  as  Adding  the  homogeneous 

part,  the  solution  in  the  inner  limit,  which  will  permit  matching  with 


t he  i nner  so  I u t i on , i s 


111  ^(1)  [ JL'  uKj  + L±  ] «. 

K*  lb  IW  L |^ /*  191  If 'I* -I 

- 0.  [ ( U lfl)‘  -—,]■> 

/</4  i(>  it*  *y  L if’i* J 

4 RP-(f)  4 


(3.17) 


This  can  be  expressed  in  terms  of  the  Inner  variables  x , * , and  y , 
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1 


cf.  (3.8b),  to  yield  finally 


* "&(  WWi 4 

*efj^  (ft  V)  * * ^ 6®  ROJ  -jff.  ♦ 

. *@  n is  fulfil  t !ki‘i  — t 

TiT;T'7f  L fe.?l  i?l’Jlfr 

* T.  [ f^D'  * **.(&)  * JpJ 4 


/6  K”J  n*  <*y 

4 i§  L'lfcJl  * Tk(”-,‘rs< ) * fy'l 


(3.13) 

where  the  first  two  terms  proportional  to  f0  results  from  expansions  of 
the  coordinates  x ' and  y"'  , cf.  (3-8b). 

Adding  the  result  of  (3-18)  to  (3-9)  gives  <$,  + in  the  inner 

. A A v 

limit.  Comparing  this  wi  th  ( <£()  f rom  (2.28)  and  (2.29),  indicates 
that  except  for  r,(y)  which  does  not  affect  the  pressure  and  the  <p£ 

^ /\  T*  ^ 

and  <p£  the  outer  expansion  of  9*  and  the  inner  expansion  of  cos ^ 
do  match,  including  a I 1 terms  exhibited  in  (2.28),  (2.29),  ( 3 • 9 ) and  (3.18), 
after  determining  n , (Tr  and  (pl  according  to  (3.10)  and  choosing  the 
doublet,  source  and  the  vortex  strengths 

EV(y)sff:rrMK.f  d:w  , 

D'  If)*  W,  (>•)  - (K.)"  Blcf)  * ««)  *-SK>-'n7  , 


Q,(y)s  £T.tr'>  = <3,00  + «*•■><>"  KSf-'ir.jfaf  >] , 
ITcrJ  = F.(?')-f,(f)  , 

(3-19) 

— A 

Except  in  the  relation  r:  (y ) =‘“Ar0g-) . ? and  y In  above^as  well 
as  in  (3.10)^are  interchangeable  without  altering  the  order  of  the  error 
^ even  for  ©»0(€_'*)  , cf.  $2J  ) . Apart  from  the  line  doublet  associated 
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with  the  lift  in  Van  Dyke's  (19&M  "third-order  theory",  the  corresponding 

solution  in  the  present  theory  has  a line  source  and  a line  doublet  associated 

— / £ -\r 

with  the  local  X -direction.  In  addition  to  the  and  D,  determined  from 

the  inner  solution,  there  is  a dominant  part  in  Q * and  p ' contributed  by 
the  nonlinear  lift  effect  proportional  to  fl*  In . 
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4.  STRAIGHT  OBLIQUE  WINGS 

For  applications  to  the  type  of  oblique  wings  considered  by  Jones 
(1971  , 1977),  the  center  line  can  be  assumed  to  be  straight,  i.e. 

© = 0,°rf®'/ « I and  PDE  (2.6)  reduces  to 

?«<  = ««©&?  (',  l) 

subject  to  errors  of  OCc*)-  If  the  variable  y is  taken  as  a normalized 
time  variable,  (14.1)  can  be  interpreted  as  one  governing  a 2-0  transonic 
sma 1 I -d i sturbance  flow  near  the  quasi-steady  limit,  in  which  the  RHS 

A A 

represents  an  unsteady  correction.  The  perturbed  solution  fo- 

under ©/  5 0 a 1 so  admi ts  * loca I -s imi lar i ty  i n ^ and  in  , appl icabl e 
to  a rather  useful  class  of  wing  geometry;  the  reduced  problems  in  this 
case  can  be  solved  once  for  a I I span  stations.  These  will  be  discussed 

i n IS  4.1  and  A. 2 below. 

4.1  An  unsteady  analogy  with  a 2-D  transonic  flow 

To  bring  out  more  precisely  the  nature  of  the  analogy  in  question, 
we  introduce 

= $ 4 «/£.  <f ‘ t , t = //«©.  (4.2) 

PDE  (4.1)  takes  the  form 

hlx-fe 3 *ftt=  . c..3) 

familiar  in  most  nonlinear  analysis  of  unsteady  transonic  flows  (Landahl 
1962,  Oswatitsch  1962,  Timman  1962).  The  equations  for  the  characteristics 
and  the  shock  relations,  (2.11)  and  (2.14),  also  go  over  correctly,  upon 
substituting  (4.2),  to  the  unsteady  transonic  system  in  question,  and  the 

A A 

continuity  requirements  on  and  , (2.10),  remain  unchanged.  The 
wing  boundary  condition  (2.9)  becomes 

(■if  Is  <*■*> 

The  far-field  description  for  If  r based  on  the  expressions 

for  £ and  £ of  (2.28)  and  (2.29),  wi th  replacing  «©%/  and @*0  's 

precisely  that  describing  the  solution  to  (4.3)  at  large|£|  near  the  quasi- 
steady limit.  One  must  note  that  the  induced-upwash  term  t.  in  if  , 
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(4.2),  cancels  out  the  corresponding  term  in  from  (2.29),  and  therefore, 
the  resultant  far-field  description  for  <P  is  independent  of  the  influence 
of  the  wake  vorticity.  Hence,  our  inner  airfoil  problem  is  mathemat- 
ically equivalent  to  that  of  a (special)  unsteady  transonic  airfoil 
pertaining  to  the  nonlinear  (lower-frequency)  domain  governed  by  PDE  (4.3) 
near  the  quas i -steady  limit,  with  a wing-surface  condition  (4.4),  and  time- 
dependent  locations  of  the  leading  and  trailing  edges.  The  wake 
influence  of  the  oblique-wing,  as  well  as  the  wing  bend  and  wing  twist, 
thus,  appear  as  additional  incidence  corrections  in  the  equivalent  unsteady 
problem. 

The  foregoing  examination  indicates  that  the  solution  to  the  oblique- 
wing problem  can  be  generated  f rom  an  uns teady  2-D  calculation  based  on  (4.3) 
for  which  existing  numerical  procedures,  e .g . those  using  alternating 
direction  implicit  algorithms  similar  to  Ballhaus  & Goorjiian  (1977),  can 
be  quite  readily  adopted.  The  crucial  input  to  such  an  application  lies, 
of  course,  in  the  incidence  correction,  and,  to  fully  utilize  the  (available) 
computer  storage  and  time  for  this  purpose,  the  description  for  J P based 
on  (2.28)  and  (2.29)  should  be  used  at  the  far  boundary.  (The  values  of 
(T0 , Dj  and  Da  in  (2.29)  can  be  taken  from  those  determined  from  the 
previous  time  step.) 

The  approach  via  this  unsteady  analogy  promises  an  alternative  solution 
procedure  without  the  requirement  for  the  local  similarity  (cf.  (4.5)  below), 
as  well  as  providing  a method  for  capturing  shock  waves  on  oblique  wings 
without  the  difficulty  associated  with  the  reexpansion  singularity  discussed 
earlier  ( § 2.6).+  Sample  calculations  via  the  unsteady  analogy  and  their 
comparison  with  the  perturbation  solutions  are  given  at  the  end  of  ^ 5 for 
a subcrltical  case.  From  the  viewpoint  of  unsteady  transonic-flow  analysis, 
it  is  of  interest  to  point  out  that  putting  D = «©»/»?  =’/*£  (2-35). 


With  the  inclusion  of  on  the  RHS  of  (4.3), the  unsteady  analogy  can 
be  used  as  a procedure  for  solving  PDE  (2.6)  with  a nonvanishing  center- 
I ine  curvature. 


the  treatment  of  the  nonuniformity  at  the  shock  root  presented  In  j 2.6 
applies  equally  well  to  the  corresponding  problem  arising  from  the  unsteady 
perturbation  considered  by  Ehler  (197**).  Traci,  Albano  & Farr  (1976), 

Hafez,  Rizk  6 Murman  ( 1 977),  Ba I 1 haus s Goorj i an  (1978),  Fung,  Yu  & Seebass 
(1978),  also  Nixon  (1978). * 

A. 2 Local  s i m i I a r i t i e s 

Owing  to  the  linearity,  the  3'0  correction  to  ^ , t $,  , can  be  de- 
composed into  suitably  scaled  separate  parts.  For  a certain  class  of  oblique- 
wing geometry,  these  separate  parts  have  similarity  solutions  independent 
of  y , as  does  the  basic  solution  90  • This  wing  class  requires,  in  addition 
to  0'iO,  ora  small©*,  that  the  basic  wing  section  at  each  span  station 
be  generated  from  a single  airfoil  shape  of  the  same  thickness  at  a fixed 
incidence.  More  specifically,  it  requires  a form  of  wing  coordinates 

♦/’Z.W'-ff  fr?>]  <s,5) 

where  C(f)  is  the  ratio  of  the  local  wing  chord  £fp)to  the  root  chord  Ca  . 
Implicit  is  that  j?rO  is  a straight  axis  on  the  wing;  the  planform,  as  well 

A 

as  the  functions  z.<?>  and  I (?)  • are  otherwise  quite  arbitrary. 

For  this  case,  we  introduce  the  variables  ( with  t = sec/\  ) 

7m  */£  *l*/c(?)  , * = , 7 CeiA  ? = 7'/t,',(l*.b) 

and  assume 

$/e=  $.  ><[[*.£ - 

<■  ••• 

where  2>  , 7.  and  A,  are  independent  of  y , and  C tdc/jy.  The  separation 

0 1 ^ 1 

into  and  <Pf  allows  the  i nduced-upwash  effect  to  be  treated  independently 
of  the  spanwi se-compress i bi 1 i ty  correction.  The  POE  (2.20)  and  (2.22) 


In  the  last  two  works  cited,  the  reexpansion  singularity  and  the  resulting 
anomaly  does  not  appear.  Its  absence  is  presumably  a consequence  of  the 
assumption  of  a pi ane , norma  1 shock  stipulated  therein. 
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become  in  this  case 


Jf  [K*  $0?  - *£  (4>ot ) ] * fjx  ft  = 0 , 

(r*.  - ««>*..-]£.  * ,4.  -«-){•.«, 4}^= - 

- 2 (*  ♦ Iff  ) & s-  ; 

I fK"  ' <<")  - ,4.  • }*<  = O. 


(4.8a) 


(4.8b) 


(4.8c) 


The  right-hand  member  of  (4  . Sts ) arises  from  9/gf  - l/py*  - (?'/ 1 ) ( tty?  + ?*/ % ) 
The  conditions  on  wing  surfaces  and  the  wake,  (2.21)  and  (2.23)  . give 


( 4 . lhi , h ,c  ) 


lt.ri-lM‘0  l%e-MttUmo  (4. io.».b.C) 

Tv  fe  re  re  rr  re 


Consistent  with  the  similarity  structure  of  $/£  . (4.7) , the  equa- 
tion of  the  shock  boundary  (2.25)  take,  the  form 

♦ *©  C'(?)Z?(C)  +*[jKK£t(f)  Z'9  l4n' 

The  Rankine-Hugoniot  relations  transferred  to  the  unperturbed  shock 

boundary,  i.e.  (2.26)  and  (2.2?) , now  read  as  at  s'  r x"P(t).+ 

o 

K„-  «*)(£*)=  ( 4 • 12a) 

= - 2p*%*)(#c%0,  (4.12b) 

- (*+0<^2-  * X,0  (4.12c) 


These  jump  conditions  for  0,  and  are  not  the  same  as  those  for  the 
(formal)  weak  solutions  of  PHE  (4.8b)  and  (4.8c). 


and 


n K i = a t ♦ a ft  ♦ * o . 


(*».  13a.b.c) 


Approaching  the  outer  limit,  I*  | =/?/*/-•>■*  • PDE  (4.8a,b)  admits 
the  developments 

% ~ $[***' (jOt) *f*i9"*l  * * fyjK*)/}?it+ 

’tz.  tK's‘/i{r]  vfci'*-  <*■«*•) 

w- 

+ ,%  *-ifl  ♦ 5 * f *»"*■  3 * 

?,  ~ ^ rta" Tc?)  * f 5»" «]  * 

♦ C'-if  i * k. 

♦ (&'.*  + &i‘)  / If  I*  * '■  • (b.uc 

These  results  are  recoverable  (completely)  from  the  far-field  expansion 

for  f t-  € £,♦•••  in  the  more  general  case,  (2.28)  and  (2.29), 

through  t rans forma t i ons  (4.6)  and  (4.7)  and  the  relations  for  the  coefficients 

JTl=K/c,  D*.“  $o/£l  , = jf)1"0)  (4.15) 

Jr' i _ Jr  * @ / ITL  ^ 1 A 

<r.  - e ■-&(?#•) c inc  • i'1-16-’1 


»©2,]r. 


(4.16b) 

(4.16c) 
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We  note  that,  according  to  the  last  of  (4.15), 


The  1 i n(2»*c  ♦ ) above  explains  the  formal  difference  of  ( hi? I * x ) i n 

(2.30)  and  '"iff  in  the  second  term  inside  the  curly  bracket  in  (4.14b). 
The  source  strength  C can  be  explicitly  evaluated  in  terms  of 

in  this  case  wi  th  the  help  of  the  Green  theorem,  unaffected  by  the 


and 


presence  of  Imbedded  shocks  (Appendix  III), 


q = *4  d:  ♦ 


j ’±1 


8"K 


(4.17) 


where  S'-  ^ and  are  the  leading  and  trai I inq-edge  locations 

in  x , respec t i ve I y . 

«N.  * 

With  £ | known,  ways  to  control  t ho  Induced  upwash  by  wlnq  twist  and 
wing  bend  through  (4.7)  are  evident.  From  the  form  of  solution  (4.7)(  the 
similitude  of  the  }-D  flow  structure  is  also  apparent . In  flow  regions 
removed  from  shocks,  the  local  pressure  coefficient  Cp  S l U* 

can  he  correlated  as 


(V-<Va-D  )/( © = f (f  K*)( 


( 4 . 18a) 


independent  of  y , where  p is  the  local  value  of  Cp'  computed  accord- 

ing to  the  2-D  theory  after  an  incidence  correction 

(4.18b) 

The  3'D  effects  on  the  field  near  the  shock  and  the  shock  boundary  itself 
may  also  be  correlated  through  (4.11)  - (4.13),  excluding  the  vicinity  of 
the  shock  root.  In  the  latter,  the  pressure  field  should  be  correlated 
in  the  original  form  of  the  Oswatitsch  f.  Zierep  (I860),  shown  in  (2.39). 
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5.  COMPUTATION  METHODS  AND  EXAMPLE:  STRAIGHT  OBLIQUE  WINGS 


5-1  Remarks  on  examples  and  comparison 

Demonstration  of  solutions  to  the  reduced  problems  of  the  theory  is 
an  essential  part  of  the  present  study.  This  is  so,  not  only  because  it 
affords  an  opportunity  for  assessing  the  adequacy  of  the  present  approach 
as  an  aerodynamic-analysis  method,  but  because  the  existence  and  uniqueness 
of  the  solutions  to  the  reduced  problems  have  not  been  thoroughly  investi- 
gated. This  demonstration  is  given  below  for  the  class  of  straight  oblique 
wings,  for  which  similarity  solutions  exist  , thus  the  sweep  and  other  3"D 
effects  can  be  more  clearly  delineated  with  the  least  geometrical  complexity. 

In  the  following,  the  computation  procedure  is  described  ( § 5-2),  shock- 
free  examples  (including  two  with  slightly  supercritical  component  flows) 
are  presented  ( § 5.3).  Comparisons  with  corresponding  solutions  from 
3-D,  full  potential  computer  programs  (the  Jameson  FLO  22)  will  be  discussed 
in  § 5.^.  A comparison  with  solutions  obtained  via  the  unsteady  analogy 
is  also  shown  in  §5-5.  Numerical  solutions  involving  imbedded  shocks, 
which  require  treatment  of  the  logarithmic  singularity  at  the  shock  root, 
are  the  subject  of  a separate  paper. 

The  present  theory  is  limited  by  the  assumptions  of  a high  aspect  ratio 
and  of  the  small  disturbance;  it  is  uncertain  without  a demonstration  that 
the  solutions  may  adequately  predict  the  (inviscid)  aerodynamic  characteristics 
to  a degree  enjoyed  by  the  classical  lifting-line  theory.  In  addition,  the 
possibility  of  committing  algebraic  errors  is  not  low.  This,  together  with 
the  two  limitations  men t ioned , makes  d i rect  comparison  with  more  exact  (3~D 
full  potential)  solutions  essential. 

The  limitation  on  the  computer  storage  available  to  the  current  3_D 
flow-field  computation  programs  is  well  known.  With  this  limitation,  it 
is  not  clear  whether  the  grid  distributions  in  these  programs  are  suffi- 
ciently refined  for  the  purpose  of  describing  the  induced  up-wash  of  the 
trailing  vorticities  of  the  far  wake,  so  crucial  to  a high  aspect-ratio 
wing.  This  uncertainty  is  erased,  however,  by  the  consistently  good 


H 


! 


i 
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comparisons  shown  below. 

5.2  Computation  procedures 

L ine-relaxat ion  methods  are  used  to  solve  the  algebraic  systems  of 
the  difference  equations  for  the  similarity  solutions  ^ , and  ; 

the  methods  employ  Murman's  type-sensitive  difference  operators,  including 
a shock-point  operator  (Murman  & Cole  1971,  Murman  197*0-  The  basic  com- 

/V* 

puter  code  for  ^ is  adopted  from  one  used  earilier  in  Hafez  & Cheng's 
(1977a  & b)  studies,  improved  in  the  far-field  description  with  (A.lAa) 
and  in  the  use  of  a higher-order  convergence  acceleration  scheme  (Meng 
& Cheng  1977).  The  double1  strength  Q*  and  in  ( 4 . 1 4a ) are  determined 

by  a least-square  fit  of  daia  near  and  at  the  far  boundary.  The  inclusion 
of  the  doublets  and  the  nonlinear  terms  shown  in  (i(.1*4a)  have  improved  sub- 
stantially the  accuracy  and  internal  consistency  of  the  numerical  solutions 
at  the  far  boundary;  this  results  in  a two  to  five  percent  (2-51.)  change 
in  the  surface  pressure.  The  basic  program  will  also  furnish  input  data 
for  the  relaxation  solutions  involving  shock- fi t t i ng  algorithms  (Hafez  & 

Cheng  1977b),  the  latter  is  not  needed,  however,  in  the  shock-free  examples 
considered  below. 

The  procedures  for  and  solving  the  linear  PDE's,  with  conditions 
on  the  x-axis  and  the  far-field  descr  i pt  ion  ( 4 . 8)  - (A . 1 6)  , are  similar  to 
that  for  . The  programs  are  simpler  in  that  the  transition  boundary 

O' 

is  fixed  to  the  sonic  boundary  of  Pe  , but  it  requires  a larger  storage 

for  the  inclusion  of  in  addition  to  the  $ S obtained  in  previous  iterations. 

0 ’ o o O' 

The  same  grid  with  nonuniform  mesh  is  employed  for  9>0  , <pt  , and  <pj  , cover- 
ing a region  lx  I £ 5 ,/C/  £ 6,  with  a total  of  81  x 65  grid  points.  The 
leading  edge  is  made  to  locate  at  x”s  -I,  and  the  trailing  edge  at  x - | 

( i . e . fl’s  - 1 , 5*  = f ) , any  departure  of  from  unity  is  accounted  for  by 


Unpublished  studies  by  Ronald  C.  Smith  at  NASA  Ames  Research  Center  appears 
to  have  raised  doubts  on  certain  earlier  results  from  the  3"D,  full  po- 
tential program  (FLO  22)  used  for  oblique-wing  analyses.  We  also  point 
out  that  currently  available  3-D  computer  programs  based  on  the  transonic 
smal l-disturbance  equation  (e.g.  Bailey  6 Ballhaus  1972)  are  inapplicable 
to  wings  without  bi-lateral  symmetry.  Therefore,  comparisons  with  3-D 
smal l-disturbance  computer  solutions  have  not  been  made. 


!i 

I 


r 
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changing  x into  [x'  + O - a ) ] in  (^4.8)  through  (^.16)  . The  case  1*1  = i 
corresponds  to  an  oblique-wing  planform  with  fore-and-aft  symmetry  about 
the  (straight)  center  line.  The  Kutta  condition  is  implicit  in  the  programs 

<V  <V  A* 

in  requiring  the  continuity  of  <P0  f ,w,  Z and  at  the  trail  i ng  edge. 

The  iterative  solutions  use  a relaxation  factor  1.3  in  the  elliptic 
region  and  0.8  in  the  hyperbolic  region;  typ  i ca 1 1 y,  200  line-sweeps  are 
needed  for  the  convergence  of  the  circulations  to  within  I0~*  • The  cal- 

A/ 

culation  for  <pt  involves  more  work,  and  requires  twenty-five  (25)  minutes 
on  an  IBM  370/158  Version  III,  using  double-precision  a r i thmet i cs .+ 

5.3  Examples  of  solutions 

As  examples,  we  apply  the  solution  procedure  to  oblique  wings  with 
the  section  function  2 *(?) genera  ted  from  the  NASA  3612-02,40  airfoil, 
scaled  to  an  arbitrary  thickness.  This  airfoil  has  been  used  in  various 
wind-tunnel  and  preliminary-design  studies  of  oblique  wing  at  Mach  number 
between  0.60  and  1.4  (Black,  Beamish  6 Alexander,  1975;  Jones  1977).  for 
the  present  purpose,  we  assume  that  the  straight  axis  (p's  0)  to  coincide 
with  the  mid  chord  of  each  wing  section;  the  spanwise  distribution  of  the 
wing  chord,  is,  however,  left  arbitrary.  In  this  appl icat io^  i t suffices 
to  set  the  ratio  of  the  thickness  and  camber  parameters  V*  equa 1 to  un i ty  , 

/V  A/ 

replacing  all  ot  by  t in  scales  entering  the  definitions  of  Z and  <p  . 
as  well  as  f(n  , @ and  £ . The  determination  of  , 9,  t and  requ  i res 
the  specification  of  the  component  similarity  parameter  K„- 

Three  sets  of  results  for  wing-surface  distribution  of  («)*.(£)!•  ■ 
and(<$J ) Z are  presented  in  Figs.  3 and  4 . The  first  set  (Fig.  3)  is  com- 
puted  for  K„-3. 6 with  Z CO  taken  from  the  rescaled  coordinates  of 
NASA  3812-02,  40  at  zero  incidence.  The  second  and  the  third  sets 

A computer  program  similar  to  that  for  Pi  has  been  considered  in  the 
context  of  a straight  unyawed  wing  by  Small  (1978). 

f-  ^ 

The  coordinates  of  the  airfoil  section  NASA  3812-02,40  may  be  expressed  as, 
with  *'s  (K'-*U>/cCr) , 

*'/lorc(y)  = (0.076/1.10')  Z’( /- *'*)  t ( o.roft'  t $'[-  o.jogr  t 

+ otzst,  * O370IK  ')}jj 
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(Figs,  4a,  b & c)  are  computed  for  a lower  similarity  parameter,  K„—  3.45, 
which  gives  a slightly  supercritical  2-D  component  flow  in  ^ ; the  airfoil 
section  and  incidence  in  the  second  set  (solid  curves)  is  the  same  as  in 
the  first;  the  third  set  (dash  curves)  differs  from  the  second  only  because 
a nonzero  incidence  4o(  = 0.0582 T (Yadian)  is  added  to  </*w/rfx',T  From  these 
results,  the  pressure  coefficient  C'S(r>-ha)/tfi,  It*  on  the  wing  surfaces 
is  calculated  from 

Cp=  coSA  Cp  z -l  7 (?,)z  + 


+ *0  ^ (%)-  t €[/iC  $f(y)  f I ♦ ®ZB]($)~j  (5.1) 

In  Fig.  3i  the  computed  surface  distribution  of  <$)r  , and 

fob  are  presented  as  thin  solid  curves,  heavy  solid  curves,  and  dash 
curves,  respectively,  with  the  itallic  and  "l  ” referr i ng  to  the  upper 
and  the  lower  surfaces.  The  peak  of(£)*  occurring  near  the  mid-chord 
point  lies  slightly  below  the  cri  t ical  va  1 ue  Kn/(4*i)~  1 • 50  . The  surface 
values  of  (£)*  which  arises  from  the  spanwise  variation  of  the  compress- 
ibility correction  are  seen  to  be  numerically  small  as  compared  to  (£)*• 
The  latter  accounts  for  the  incidence  correction  and  includes  the  far-wake 
vorticity  influence.  The  circulations  given  by  the  jump  of  P0  , , and 

^ at  the  trailing  edge  are  found  to  be  Jf^  — 2 . 0^3 , /]T  r -0 . b3b , and 
— -2.667,  respectively.  Their  contributions  to  the  potential  9 at 
the  trailing  edge,  hence,  the  rolling  moment  about  the  wind  axis,  are 
weighted  by  C > €®  C , and  ]>  respectively.  For  a 

pi  anform  which  is  symmetrical  with  respect  to  y'  — 0 (not  y=0)tff'o  is 
also  symmetric.  The  induced  upwash  in  this  case  is  positive  C <0)  on 
an  aft  panel  and  negat  i ve  ( (E,  >0)  on  a fore  panel,  as  noted  earlier. 
Therefore,  the  last  two  of  the  three  weighting  factors  mentioned  are  both 
negative  on  an  aft  panel  and  positive  on  a fore  panel.  It  is,  thus,  seen 
from  this  example  that,  with  the  negative  and  [J^  , the  inviscid  3-D 
effects  will  give  rise  to  an  unbalanced  rolling  moment  (as  well  as  pitching 

T For  X - 0 . / 2 , this  amount  of  A oC  gives  an  incidence  of  0.40  degree. 
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moment).  The  asymmetrical  forces,  if  unchecked,  will  tend  to  raise  the 
aft  panel  and  lower  the  fore  panel. 

The  least-square  fit  of  the  doublet  strengths  determined  from  the 
relaxation  solution  of  Fig.  3 are  &r0  = 0.2972  and  — 0.0958; 

with  these  and  the  value  of  IT0  — 2.098,  the  source  strength  in  the  far 
field  for  <r(  is  determined  from  (9.17)  as  = 0.9270.  We  note  in 

/v 

passing  that  the  program  for  computing  is  relatively  straightforward 

/V 

and  the  result  can  be  reproduced  quite  well  by  computing  4>0  at  a slightly 
different  incidence,  taking  the  difference  of  the  two  j , and  normalizing. 

To  assess  the  relative  importance  of  those  contributions,  we  apply 
the  results  of  Fig.  3 to  an  elliptic  planform  with  a major- to-m i nor  axes 
ratio  of  16.78  at  30°  yaw.  Assuming  a 67.  thickness  ratio  for  each  wing 
section,  we  have  for  this  example  £ = 0.1522,(5  - 1.337  (corresponding  to 
Mb=  0.7615  , and  M.=  0.8793  under  Kn~  3-6).  No  wing-bend  is  considered, 
but  a uniform  If?)  corresponding  to  an  (overall)  incidence  adjustment  is 
chosen  to  eliminate  3'0  effects  on  the  total  lift  (this  is  possible  owing 
to  the  pure  antisymmetry  of  the  3'D  effects  in  this  special  example).  The 

A. 

resultant  distribution  of  0,1  the  upper  surface  at  the  span  station 
y — 0.80  is  shown  as  a dotted  curve  in  Fig.  3.  A small  region  of  super- 
critical component  flow  is  seen,  which  in  this  case  is  brought  about  mainly 
through  the  upwash  induced  by  the  wake  vorticities. 

The  corresponding  results  for  the  other  two  sets  shown  in  Figs.  9a,  b 
and  1 are  similar,  with  the  noticeable  difference  being  the  appearance  of 
supercritical  regions  on  the  upper  surface  in  the  leading  approximat ion 
(«)?.  Tlie  computation  does  not  experience  convergence  difficulty  and 
examination  of  the  solution  does  not  reveal  evidence  of  shock  waves  in 
these  cases.  For  the  second  set  with  Kn  = 3-95  anil  A oi  - 0 (zero  incidence), 
the  solutions  give  =•  2.0712,  = -0.5291  , fl\=  -2.6839,  dJ  = 0.2599, 

D*’  — 0.0071  and  Q = 0.9530;  for  the  third  set  with  1C  =r  3-95  and  an 
incidence  of  0.9  corresponding  to  A*-/x  = 0 .0582,  the  cor respond i ng  results 
are  = 2 . 2828 , |f’  = -0.5927,  = -2.2118,  = 0.2619,  D‘  = 0.02578. 

and  = 0.9893. 


5 . 5 Compa r 1 son  w 1 1 h fu I I -potent i al , 3~D  sol ut ions 

2^  ^ a. 

The  similarity  solutions  , t,2  and  shown  in  Figs.  3 and  5 can 

be  applied  to  examples  of  transonic  oblique  wings  of  various  combinations 
of  \^  , T , A as  long  as  the  component  similarity  parameter  K„  is  unchanged. 
The  wing  aspect  ratio  and  the  spanwise  distribution  of  the  local  wing  chord 
can  be  quite  arbitrary. 

Surface  distributions  of  Cp  have  been  obtained  in  several  cases 
for  the  purpose  of  comparison  with  solutions  from  a computer  program  based 
on  a 3-D  fu 1 I -pot en t i a I equation.  In  this  connection  it  is  essential  to 
differentiate  the  significance  of  the  critical  pressure  coefficient  corre- 
sponding to  a local  sonic  speed, 


and  that  corresponding  to  the  sonic  condition  in  the  component  flow,  namely. 


(5.3) 


For  the  viewpoint  of  swept-wing  aerodynamics,  Cp  is  relatively  trivial, 
inasmuch  as  a three-dimensional  flow  with  Cp  < Cp*  can  readily  be  rendered 
shock-free  by  providing  a sufficiently  large  wing  sweep.  The  value  of 
Cp4*  of  (5.3)  is  very  well  approximated  by  that  deduced  from  the  present 
study 

The  difference  between  (5.5)  and  (5.3)  apparently  depends  on  lzJ  ^ ^ , 

which  are  found  to  be  less  than  2 for  as  low  as  0.65. 

It  may  be  pointed  out  that  existing  3~D  computer  codes  based  on  the 
ful I -potent ial  equations  do  not  account  (correctly)  for  the  departure  of 
the  trailing  vortex  sheet  from  a planar  (flat)  surface.  Thus,  its  validity 
also  requires  the  sma 1 1 -d i s t urbance  assumption,  strictly  speaking.  The 
value  of  the  program  lies,  of  course,  on  its  ability  to  describe  the  flow 
field  around  the  leading  edge  and  near  the  trailing  edge,  where  the  sma I I - 
disturbance  approx  I mat  ion  breakdown  or  becomes  less  accurate,  assuming 
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that  problems  of  convergence  with  respect  to  mesh  size,  and  to  iteration, 
may  not  cause  inaccuracy. 

Numerical  results  comparable  with  our  solutions  are  generated  from 
one  version  of  A.  Jameson's  3"D  fu) 1 -potent ial  computer  codes  "FLO  22" 

(c.f.  Bauer,  Garabedian,  Jameson  6 Korn  1974,  Jameson  1974),  which  is  used 
with  implementations  for  oblique-wing  analyses  at  NASA  Ames  Research  Center 
Ae ronaut i cs D i v i s i on  and  at  Grumman  Aerospace  Corporation  Research  Depart- 
ment. The  algorithms  employed  in  FLO  22  are  not  fully  conser va t i ve , but 
this  may  not  be  essential  for  shock-free  solutions  presented  below.  We 
point  out  that  the  FLO  22  data  from  Ames  and  from  Grumman  are  not  identical 
owing  mainly  to  the  use  of  different  meshes.  The  availability  of  data 
from  two  sources  is  helpful  in  delineating  the  nature  of  discrepancy  between 
our  theory  and  the  more  exact  3"D  program.  Data  from  t fie  latter  is  still 
influenced  by  the  mesh  size,  spacing  of  the  span  stations,  number  of  iterations, 
the  detail  of  the  leading-edge  geometry  description,  which  are  different  in 

i 

Ames  and  Grumman  runs. 

A number  of  FLO  22  runs  have  been  made  with  free-stream  Mach  number, 
swept  angle,  wing-thickness, etc . chosen  to  give  either  K„=  3-6  or  — 3.49, 
employing  the  same  basic  airfoil  section.  An  elliptic  planform  is  used  in 
each  case;  wing  twist  and  wing  bend  are  assumed  to  be  zero.  Comparisons  are 
made  in  Fig.  9“7  for  three  cases  and  may  be  considered  being  typical  among 
most  runs. 

The  oblique  elliptic  wing  considered  in  Fig.  5 has  a 6‘7.  thickness  ratio, 
with  the  major- to-minor  axes  ratio  of  20,  and  a sweep  angle  of  22.5°.  The 
free-stream  Mach  number  is  0.0242  in  this  case.  This  makes  the  free- 

stream  component  Mach  number  = 0.7619  and  K„=  3-60, (0  = 1.003  and 

£ = 0.1277  in  our  analysis.  The  pressure  coefficient  for  the  critical  com- 
ponent  flow  in  this  case  is  Cp  = -0.470.  The  surface  Cp  values  at  three 
span  stations,  y =-0.69,  — 0 and  — 0.69  are  shown  in  Fig.  9a,  Fig.  9b 

and  Fig.  9c,  respectively.  The  FLO  22  data  from  Ames  (in  small  crosses  and 
'•y  " and  from  Grumman  (in  small  open  circles)  are  seen  to  be  quite  close 


+ The  FLO  22  data  generated  at  the  NASA  Ames  Research  Center  were  kindly 
provided  by  Mr.  Ronald  C.  Smith,  and  the  data  from  Grumman  Aerospace 
Corporation  were  made  available  to  us  by  Dr.  Rueben  Chow. 
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except  near  the  leading  edge.  Near  the  leading  edge,  the  sma I I -d i s t urbance 
assumption  of  our  theory  al so  breaks  down  as  noted  earlier.  The  agreement 
of  the  FLO  22  data  with  our  values  computed  on  the  basis  of  data  from  Fig.  3 
(in  solid  curves)  must  be  considered  as  being  better  than  anticipated, 
inasmuch  as  the  relative  error  in  our  theory  is  of  an  order  determined  by 
the  larger  of  and  €*  (to  be  precise).  In  the  present  case,  X ^ = (0.06) 

= 0.153! 

Encouraging  is  that  the  degree  of  agreement  with  the  FLO  22  data  appear 
to  deteriorate  little  with  increasing  wing-thickness  ratio  or  reducing  aspect 
ratio.  Results  of  comparisons  for  elliptic  wings  with  a 12V.  thick  airfoil 
section  and  a major- to-mi nor  axes  ratio  of  1 4 are  shown  in  Figs.  6 and  7. 

The  case  considered  in  Fig.  6 has  a free-stream  Mach  number  0.7549,  a yaw 
angle  30°.  This  gives  0.6538,  K„  = 3.60,  6 s 0.1448,  and  (*)  s | .062; 

a I so  cr = -0.727.  The  Cp  values  (compted  from  the  present  theory,  based 
again  on  data  of  Fig.  3)  are  presented  for  five  span  stations  in  heavy  curves 
with  open  circles;  the  corresponding  FLO  22  data  (from  NASA  Ames)  are  pre- 
sented as  thin  curves.  Agreement  similar  to  that  found  in  Fig.  5 is  again 
found  in  Fig.  6 except  at  the  station)  y = £0.89  which  are  rather  close  to 
t he  w i ng  tip. 

In  Fig.  7.  the  thickness  ratios  and  aspect  ratio  remain  the  same  as  in 
Fig.  6,  but  0.7677  and  A*  30°.  This  gives  M - 0.6648,  Ce**  - -0.689, 

and  K = 3.3*5.  For  this,  the  C«  value  in  the  present  theory  is  computed  on 
the  basis  of  our  data  from  Fig.  4.  The  results  are  presented  in  Fig.  7 for 
seven  span  stations;  y * 0,  1 0.20,  i 0.59,  i 0.79.  The  degree  of  agree- 
ment with  the  FLO  22  data  (from  NASA  Ames)  is  again  similar  in  the  preceding 
comparisons.  In  this  case,  a small  supercritical  region  appear  in  the  aft 
wing  panel,  as  anticipated  (there  are  at  least  three  surface  grid  points  in 
the  supercritical  region). 


+ The  two  sets  of  FLO  22  data  also  differ  near  the  wing  tips  ( y~  z * OlO) 
not  shown. 

-f  -f 

Note  that,  there  are  considerably  more  data  points  than  the  numbers  of 
the  open  circles  shown. 


The  reader  is  called  to  the  attention  that  the  Cp  d i st r i but  ions  are  not 
the  same  in  different  span  stations,  therefore,  the  comparisons  made  with 
the  FLO  22  data  have  provided  a meaningful  test  for  the  theory.  One  of 
the  major  contributions  to  the  3*0  effects  is  the  incidence  correction  which 
accounts  for  the  upwash  induced  by  the  wake  vorticities,  computed  from 
Y*of  (3-9)  via  (3-10)  and  (4.16a).  The  adequacy  of  this  upwash  calcu- 
lation has  already  been  demonstrated  in  examples  pertaining  to  the  linearized 
problem  (Cheng  1978a  £ b) ; therefore,  the  good  agreement  reported  here  should 
not  be  too  surprising.  Worthy  of  note  in  this  connection  is  the  use  of  the 
transonic  similarity  parameter  K„=  0 , and  the  retention 

of  the  cos_/\  factor  in  computing  and  Cp  , which  prove  to  be  crucial  in 

maintaining  an  accurate  leading  approximation,  and  is  partly  responsible 
for  the  good  agreement  achieved  here. 

5.9  Comparison  with  sol ut ions  via  the  unsteady  ana  I ogy 

To  demonstrate  the  degree  to  which  the  unsteady  transonic  analogy  can 
be  utilized  for  the  oblique-wing  analysis,  we  include  here  a comparison  of 
the  local -s imi lar i ty  solutions  for  non- 1 i f t i no  oblique  wings  with  the 
corresponding  results  generated  by  a numerical  procedure  using  an  ADI  algo- 
rithms similar  to  those  in  Ballhaus  6 Goorjian  (1977).+  As  pointed  out  in 
§ 4.1,  the  case  involving  lift  ( (T0  £ Oi)  regu  i res  an  implementation  at  the 
far  boundary  for  the  unsteady  analogy,  not  needed  for  the  example  considered 
here.  The  result  is  nevertheless  of  interest  in  that  the  3*D  correction, 
in  this  case;is  given  entirely  by  the  compressibility  effect  resulting  from 
spanwise  density  variation  mentioned  earlier.  For  this  purpose,  we  consider 
oblique  wings  with  a symmetric  airfoil  section^NACA  64006  at  zero  incidence, 
and  compute  ^ for  = 2.50;  (cf.  § 5.2)  the  result  for  the  surface  pressure 
coefficient  Cf'(#)=  was  computed  for  T = 0.06  and  0.827 


The  unsteady  transonic  code  used  in  this  study  was  made  by  Or.  Tom  Evans 
during  1977  on  leave  from  the  University  of  East  Anglia,  Norwich, 

England . 
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and  is  shown  in  Fig.  8 as  a heavy  solid  curve. + The  basic  component  flow 
in  this  case  has  a -0. 362*4,  very  closely  approached  by  the  Cp  minimum 

shown.  The  solutions  to  the  reduced  problem  of  the  3_D  correction  have  been 
computed  for  three  cases  (cases  A,  0 and  C)  corresponding  to  three  different 


chord-wise  locations  of  the  straight  axes  (center  line);  the  results  for 
the  (f.)t  on  the  surface  are  shown  as  thin  solid  curves  in  Fig.  8 (the 
maximum  value  of  in  case  C is  finite  but  exceeds  the  boundary  of  the 

figure).  The  noticeable  differences  among  cases  A,  B and  C,  confirms  the 
expectation  that  increasing  the  (local)  leading-edge  swept  angle  reduces 
the  maximum  local  Mach  number, and  decreasing  the  sweep  increases  it.  This 
is  because  the  contribution  of  $>,  X to  is  weighted  by  e © dc/ dy~  y 

the  latter  is  positive  on  the  fore  panel  and  negative  on  the  aft  panel. 

Thus  moving  the  straight  axis  to  the  leading  edge  (say,  changing  from  case 
A to  case  B) , the  $>x  maximum  will  increase  in  the  aft  panel  and  decrease 
in  the  fore  panel,  whereas,  moving  the  straight  axis  to  the  trailing  edge 
(say,  changing  from  case  A to  case  B),  the  <p  J peak  will  decrease  in  the 
aft  panel  and  will  increase  in  the  fore  panel. 

To  test  the  unsteady  analogy  in  these  (nonlifting)  cases,  we  consider 
high  aspect-ratio  obligue  wings  with  the  same  component  similarity  parameter 
( K„ = 2.50),  choosing  an  aspect  ratio  and  sweep  angle  to  correspond  to 
€@  = 0.10.  Surface  pressure  coefficients  Cp  are  computed  from  fiojf  l1nd 
$iZ  ° ^ Fig.  8 for  four  different  planforms  (an  ellipse,  a lenticular  shape, 
and  two  ha  I f- lent icular  shapes),  assuming  "t  = 0.06  and  = 0.827  for 
convenience.  Fig.  9 presents  Cp  vs.  x's(x'-ct)/2  on  a yawed  elliptic 
planform  at  three  representative  span  stations  y = -0.50,  0,  and  0.50. 

The  data  computed  from  the  present  theory  (in  dots,  crosses  and  triangles) 
compare  well  with  the  corresponding  data  computed  from  the  unsteady  transonic 


+ The  more  general  result  for  f independent  of  T and  Vl^  may  neverthe- 

less be  recovered  by  interpreting  Cp ^ in  Fig.  8 as  Cp^/z 


++ 


If  T z0 . 06  and  A - 22 . 5°,  these  values  of  K„=  2.50  and  C0  = 0.10  will  mean 
= 0.827,  (VV»=  0.895  and  ^,=  22.57;  if  X =0.12, A = 22. 5°.  then 
= 2.50  and  s 0.10  require  M^*0. 7^1.  = 0.8025  and  = 1 6 . 1 *4 . 

The  factor  0.909  in  the  ordinates  of  Figs.  9 and  10  results  from  a peculiar 
factor  introduced  in  the  Cp  of  T.  Evans'  program. 


r 
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code  of  T.  Evans  (in  thin  solid  curves).  In  Evans'  work,  the  variables 

A,  A / A # A f i /)  , 

t , * , r , Z and  r with 

t'-[l  * t/ f 0j  /l-(l4-r  )/z  ; *'=(*- °VZ,  Ts*/Z) 

are  used.  Discrepancies  are  noticeable  near  the  leading  edge,  but  are 
attributable  to  the  cruder  mesh  near  the  leading  edge  used  in  computing 

/V  A 

4>0  fr  and  <p,Z  of  Fig.  8. 

Similar  comparison  for  a yawed  lenticular  pi  an  form  are  given  in 
Fig.  10a  for  the  aft  panel  at  span  stations  'y' =.  0,  0.5,  and  0.8,  and  for 
the  fore  panel  at  span  station  y =-0.20  and  -0.50  in  Fig.  10b.  Figs.  11a 
6 b present  corresponding  results  for  a ha  I f- lent icular  planform  with  a 
straight  leading  edge;  Figs.  12a  6 b shows  results  for  the  half-lenticular 
planform  with  a straight  trailing  edge.  The  degree  of  agreement  between 
our  calculations  based  on  the  local  similarity  and  the  data  from  unsteady 
analogy  shown  in  Figs.  9 through  12  are  very  consistent.  The  encouraging 
results  appear  to  suggest  that  the  analogy  should  hold  even  for  values  of 
60  considerably  greater  than  0.10. 

Comparisons  with  unsteady  solutions  involving  circulations  and  im- 
bedded shocks,  as  well  as  the  extension  of  the  analogy  to  allow  for  center 
line  curvature,  remain  to  be  considered  in  future  works. 


6. 


CONCLUDING  REMARKS 


In  summary,  we  have  presented  a t ranson i c- f 1 ow  theory  of  thin  wings 
of  high  (reduced)  aspect  ratio  ( C « I ) which  permits  delineation  of  the 
influence  of  wing  sweep,  center-line  curvature,  and  other  3~D  effects  on 
the  nonlinear  flows  more  explicity  than  available  methods.  In  the  domain 
of  K„  and  (m)  considered,  the  outer  3'D  flow  has  a transonic  free  stream 
(Ai-f  ),  but  the  disturbance  is  weak  enough,  and  the  wing  aspect-ratio 
high  enough,  to  allow  a lifting-line  description  in  Prandt I -Glauert  vari- 
ables (as  a leading  approximation  to  ) . In  the  inner  region  next  to  the 
wing  section,  the  (local)  sweep  is  large  enough  to  avoid  strong  shock  losses 
but  low  enough  to  sustain  a transonic,  mixed  component  flow.  The  3"D 
corrections  to  this  nonlinear  component  flow  is  analyzed  by  solving  a 
perturbation  problem  for  6 1,  and  matching  the  result  to  the  outer 

solution. 

The  basic  flow  around  each  wing  section  is  the  2-D  component  flow 
characterized  by  a similarity  parameter  Kn  which  is(in  turn( control  I ed  by 
the  local  wing  sweep,  .A  . For  a given  wing  geometry,  3"D  corrections  for 
the  domain  analyzed  come  mainly  from  two  sources:  one  arises  from  a 3~D 
compressibility  correction  to  the  PDE,  absent  in  the  classical  theory,  and 
is  proportional  to  €(*)  ; the  other  corresponds  to  the  familiar  induced 
upwash  corrections,  including  the  contribution  from  the  near  wake  and  the 
velocity  induced  by  the  curved  lifting  line  itself.  The  upwash  correction, 
explicitly  determined  as  a functional  of  the  circulation  in  (3.9^)  and(3-10), 
consists  of  terms  proportional  to  € as  well  as(*)€i»i€.  On  a straight 
oblique  wing,  corrections  from  both  sources  generally  tend  to  make  the  flow 
more  critical  on  the  aft  panel  and  less  critical  on  the  fore  panel.  Never- 
theless, center-line  curvaturefas  well  as  the  local  leading-edge  swept 
angle(may  be  used  to  control  this  undesirable  3~D  effect. 

For  oblique  wings,  of  which  the  center  line  of  the  planform  can  be 
taken  as  being  straight  (®5  0 ),  the  3-0  flow  structure  possesses  simi- 
larities, permitting  the  reduced  PDE  systems  to  be  solved  once  for  al 1 
stations.  This  similitude  has  a requirement  on  the  wing  geometry,  namely, 
the  wing  ordinates  be  generated  from  the  same  airfoil  section  with  a fixed 
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thickness  ratio  and  incidence,  allowing,  however,  wing  twist  and  wing  bend 
in  the  manner  indicated  in  ( ^ . 5 ) • The  origin  of  the  wing-section  coordinates 
in  this  similitude  must  coincide  with  the  (straight)  center  line,  the  plane- 
form  is  otherwise  arbitrary.  For  straight  oblique  wings,  an  analogy  also 
exists  with  a 2-D  transonic  flow  under  the  influence  of  slow  (unsteady) 
changes  in  airfoil  geometry  and  in  the  condition  at  the  far  boundary.  The 
analogy  (demonstrated  in  the  nonlifting  cases  in  ^ 5-5)  promises  an 
alternative  method  for  oblique-wing  analysis,  utilizing  existing  numerical 
procedures  for  unsteady  2-D  transonic  flows,  without  requiring  the  local 
similarity  as  well  as  treatment  of  the  shock-root  nonuniformity  ( § 2.6). 

Such  a procedure,  with  the  upwash  function  (£*,  Cy)  determined  by  (3-9b)  and 
(3-10),  may  also  be  extended  to  the  more  general  case  with  a curved  center 
line. 

Of  interest  is  a source  term  appearing  in  the  far-field  expansion  of 
the  inner  solution,  resulting  from  the  3'D  compressibility  correction  and 
the  nonlinear  lift  effects,  this  gives  rise  to  a line  source  in  the  second 
approx i ma t i on  to  the  perturbation  potential  in  (fie  outer  flow.  It  is 
pointed  out  that,  unless  <*1  = 0,  a line  source  must  similarly  arise  in  a 
linear  theory.  The  existence  of  the  source  term  with  a nonvanishing 
strength  in  the  present  theory  is  firmly  established  for  oblique  wings  by: 

(i)  demonstrating  (in  ^ 3-3)  the  matching  of  the  inner  and  outer  solutions 
to  all  terms  larger  and  comparable  to  the  source  term  in  the  overlapping 
region;  (ii)  use  of  the  result  of  Green's  theorem  to  show  (in  Appendix  II. l) 

/V 

that  the  source  strength  Q(  , is  the  same  as  that  for  a Poisson  solution; 

( i i i ) explicit  eva I ua t ion  of  Q for  the  Po i sson  solution  (in  Append i x III). 

For  flow  involving  imbedded  shock,  the  formal  perturbation  analysis 
leads  to  a nonuniformity  at  the  shock  root,  owing  to  the  reexpansion  singu- 
larity. A nonuniform  problem  of  the  same  nature  has  been  present  also  in 
a number  of  existing  perturbation  analyses  in  the  context  of  unsteady  transonic 
flow,  as  noted.  In  the  treatment  presented  in  £ 2.6,  it  is  shown  that,  in 
coordinates  fixed  to  (and  moving  with)  the  shock  root,  the  Oswa t i t sch-Z i erep 
(I960)  description  for  the  behavior  of  pressure  or  tangential  velocity 
(not  the  velocity  potential)  remains  valid,  irrespective  of  the  perturbation. 
Thus  the  3"D(a?  well  as  the  analogous  time  dependent)  effects  can  influence 


-J 
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the  pressure  behavior  in  question  only  through  changes  in  the  three 
parameters  (local  shock  strength,  shock-root  position  and  another  constant) 
which  specify  the  Oswa t i tsch-Z i erep  singularity.  These  changes  can  be 
determined  by  matching  with  the  solution  removed  from  the  shock  root. 

The  procedures  for  solving  the  similarity  solutions  to  the  reduced 
problem  have  been  described;  solutions  are  demonstrated  for  oblique-wing 
geometries  using  a rescaled  NASA  3^12-02,^0  airfoil  section.  The  solutions 
give  maximum  flow  speeds  which  are  slightly  below  and  above  the  critical 
values.  The  results  on  surface  pressure  have  been  applied  to  specific  3'D 
oblique  wing  with  elliptical  planforms  pertaining  to  various  combinations 
of  free-stream  Mach  number,  wing  sweep,  th i ckness- ra t io  and  aspect-ratio^, 
and  compare  (in  § 5-M  with  the  corresponding  data  obtained  from  a full- 
potential,  3*0  computer  program  (FLO  22).  The  d i screpanc ies  found  are 
generally  smaller  than  the  expected  errors  from  the  sma 1 I -d i sturbance 
approximation  and  the  descret i zat ion  in  some  instances,  the  discrepancies 
could  be  attributed  to  the  convergence  problem  in  the  FLO  22  data  The 
encouraging  agreement,  cons i s tent  1 y found  in  most  cases  examined  (including 
thickness  ratios  as  large  as  0.12  and  aspect  ratios  down  to  10)  has  con- 
firmed our  expectation  from  a similarly  encouraging,  previous  comparison 
for  the  linear  theory  (Cheng  1973a  6 b) . Essential  in  this  regard  is  the 
retention  of  the  cosy\  factor  in  Jfn  for  the  inner  solution  and  in  Cp , 
and  the  adaptation  of  an  H,"  factor  in  the  similarity  parameter  K„  to 
better  approximate  the  critical  speed  ( $ 2.3).  As  yet  to  be  demonstrated 
r-  similarity  solutions  with  a more  extensive  region  of  supercritical  com- 
pcr.ent  flow,  with  or  without  an  imbedded  shock. 

For  numerical  solutions  involving  shocks,  the  type-sensitive  difference 

/V' 

schemes  are  implemented  by  shock- f i t t i ng  algorithm  for  <pg  , and  by 
corresponding  algorithm  for  ^ and  <Pi  . The  latters  are  supplemented  with 
a treatment  of  the  shock-root  nonuniformity.  Their  comparison  with  the 
corresponding  unsteady  solutions  through  the  analogy  should  furnish  an 
opportunity  to  test  the  adequacy  of  the  proposed  method  and  the  ADI  pro- 
cedure in  describing  shocks.  In  this  connection,  a solution  via  the  unsteady 
analogy  involving  lift  is  yet  to  be  carried  out.  These  works  are  in  progress. 
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In  future  work,  utilization  of  the  unsteady  analogy  should  be  explored 
further  for  its  potentiality  in  solving  swept-wing  problems  without  the 
assumptions  of  local  similarities  and  a straight  center  line.  The  re- 
laxation of  the  sma 1 I -d i s t urbance  assumption  would  appear  to  be  a natural 
extension  of  the  present  work.  It  must  be  pointed  out,  however,  that  the 
POE  governing  the  perturbation  velocity  potential  in  such  cases  is  not 
only  more  complicated,  but  does  not  readily  yield  a similarity  solution  as 
in  the  present  theory.  For  practical  consideration,  one  may,  of  course, 
replace  Qo  by  the  corresponding  ful 1 -potent ial  va 1 ue , app 1 y i ng  e<p,  from  our 
theory  as  a 3'D  correction. 

One  drawback  discourages  applications  of  the  present  work  to  conventional 
swept  wings  is  the  nonuniformity  near  the  apex  (cf.  (3 .13b)),  removal  of 
which  would  call  for  the  use  of  the  nonlinear  3'D  transonic  smal 1 -d i sturbance 
eguation  therein.  The  present  theory  may  still  provide  reasonable  results 
in  regions  removed  from  the  apex,  and  a direct  comparison  with  more  exact 
numerical  solutions  (Jameson  £.  Caughey  1977,  or  Bailey  or  Ballhaus  1972) 
for  symmetric  swept  wings  should  be  made.  The  problem  of  a supersonic  (or 
slightly  supersonic)  oblique  or  swept  wing  with  a transonic  component  flow 
represents  an  obvious,  and  equally  interesting  area  of  extension;  in  this 
case, the  problem  of  the  apex  nonuniformity  may  become  more  tractable. 
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APPENDIX  I 


THE  CRITICAL  SPEED  AND  THE  TRANSONIC  SIMILARITY  PARAMETER 

We  consider  a two~d i mens i ona I steady  potential  flow.  The  transonic  small- 
disturbance  approximation  leads  in  this  case  to  the  familiar  PDE  for  the 
perturbation  velocity  potential 

[ 0 - *£  ) - CYtOK,2  ( 1 • 1 ) 

* ys  A 

or,  in  terms  of  the  variables  x,  z and  Q defined  in  the  text  , writing 

_ . A A 

T for  £x  and  y for  c , 

f = o , (1.2) 


whe  re 


Ks, 


_ ) -M.‘ 


VC'1  TV< 


(1.3) 


is  identified  with  the  similarity  parameter  of  Sprieter  (1953).  except  for 

a change  in  the  algebraic  sign  and  the  omission  of  the  factor  C-T+/  yr/s . 

2 V/? 

As  an  asymptotic  theory,  (-/  + /)  Mi  and  /W^  in  (l.l)  and  (1.3)  may,  of 
course,  be  replaced  by  ( -f  t 1 ) and  unity,  respectively;  similarly,  ( i-*£) 
by  2(  ( -N^)  . We  also  note  that,  according  to  (2.1)  in  the  text,  (Y+l)M^ 
in  (l.l)  could  have  been  written  as  >*£  [ i - K-n)  ; but  the 

resulting  differences  are  numerically  small  in  the  range  of  considered 

( see  be  I ow) . 

The  full  potential  PDE  in  this  case  is 


(a1-  +2uv$vy+(a‘Vt)&,y=  O, 


(i.<0 


where  CL  is  the  sound  speed,  U - $ K . lr  - and  § = if  X + <P  . 
form  of  (I.M,  valid  in  the  domain  K =oco  and  x sma I I , 

l r'-by)  2 

{ I - *C-  M‘[:  * F -M’  [*? ' t 


An  approximate 
s (wi th  u'  -ft 


<Pn  0 

(1.5) 
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- ' 


This  is  one  order  X more  accurate  than  (|.l).  The  characteristics  of 
(l-M  and  (1.5)  give  the  same  critical  values  for  ( “Vtf  ) up  to  the 
second  order  of  ( 1 ~M*) 


— /v  I 

0 ~ n;  l 


I -M# 
2(H0 


(1.6) 


where  the  leading  term  ( » -m£)/(  i t 1 )*£  is  the  critical  speed  based  on 
PDE  (1.1). 

Comparing  (1.1)  with  (1.5).  the  smal 1 -di sturbance  equation  (|.l)  is 
seen  to  represent  the  characteristics  system  poorly  in  region  near  the 

critical  speed,  g i v i ng  er  roneous  s I ope  dyc/dx  therein.  The  sma  1 1 -d i sturbance 
equation  can  be  written  as 

(i-n)wj  - -jj  J 4>  ♦ <p I = 0 (1.7) 

[ «*>  V J **  • 


over  which  an  improvement  in  this  regard  can  be  gained  by  writing 

wh1^[(£)4-£K*^=  o.  »••> 

with  (U/U)#  given  by  the  two-term  expansions  (1.6).  This  modification 
then  leads  to 


[K's-  a+o$i]  $st  *■  = o , 


wi  th 


r,  _ uni  i 

frC  t:Vj  ^ -* 


(1.9) 


(l.io) 
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i 


This  expression  of  Kj  is  equivalent  to 


noting  that 


= ZV>"«  , J«H) 

k'/k;  = i»opi-nJ]. 


( i.ii) 


I 


I 

* 

i 
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appendix  I I 


APPLICATION  OF  GREEN'S  THEOREM  AND  THE  FAR-FI ELD  BEHAVIOR 

III  A relation  for  the  weak  solution 
Let  p denote  a funct  ion  satisfying 


(it'  * 


(i  i . i) 


in  an  infinite  domain  in  x and  y , , where  the  distribution  £ 

var.es  with  x and  7 and  may  also  depend  on  <p  and  its  partial 
derivatives.  Let  the  inner  boundary  of  be  the  limiting  boundary 

enclosing,  and  approaching,  the  two  segments  W (wing)  and  T (trace)  on 
the  x-axis.  also  the  segment  D (shock  discontinuity)  as  shown  in  the 
Sketch  (Fig.  A.I).  We  assume  that  across  T the  function  p satisfies 

ff  M =1TM  = o (ll.2) 

where  £ JJ  stands  for  the  jumP>£  )+  - ( )\  and  that, 

the  jump  takes  on  values  consistent  with  weak  solutions  to  (II. I), 
namely,  at  x - x°(  >)  ) , 

= o,  m = o I,,.,, 

We  assume  further  that  V * vanishes  at  the  infinity  in  the  mean,  i.e., 

- |»m  1 


across  D, 


7 p s 

R+eo  H 


1 V 


(II.*) 


where  R is  the  radius  of  a circle  centered  at  (x . fj ) and 

and  that  the  integral  of  £ £ over  ^ exist*.  Application  of  Green’s 
theorem  to  this  system  yields  a relation 


4>  (■*'?)=  j'  { Ur  - n+n  rf,  V,  Ur']d$t  + 

U/*  T f ♦ D' 


(n.  5) 


49* 


This  is  equivalent  to  (2.114),  for 


S ^ f Qy  * «©  X ' 
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/ 2 2 

where  r'  = J (x-x  J + ( fj  - rf, ) , and  O is  the  inward  normal  into 
49t.t  the  integration  is  taken  only  along  one  side  of  the  boundary 
segments.  The  last  integral  over  can  be  evaluated  as 

/i’lj.Z)  l*r  -j  TZ  Ul*r  ~ Jf  £ "~T  dx,*?< 


also. 


/ ]/*"'■'</  s,  = [iLK-tifT?  Tiff'd?,  i 

P'  ' D* 


Thus,  applying  the  first  of  (11.3), 


J. 

in 


(M.6) 


; f ( Ur'-a*l}-n  l*r'\s,  + ± l”1- 

Pr  <cK. 

We  now  apply  (11.6)  to  the  last  integral  of  (11.5)  and  the  second  of 
(11.3)  to  the  line  integral  over  D ; then,  we  integrate,  by  parts  (with 
respect  to  x,)  , the  line  integral  over  W and  T in  (11.5),  and  finally 
involk  (11.2);  we  recover  the  basic  relation  for  the  weak  solution  of 
(ll.l),  first  given  in  Murman  and  Cole's  (1971)  original  paper 

4> = ; *»- '(*?)  * #■*• 


Lt 


Lf 


- -L  ((  **Z«. £■**,</» 
in  J J r,  I 'I  , 


(11.7) 


where  the  first  two  integrals  are  taken  along  the  x-axis.  The  above 
result  still  holds  if  there  are  more  than  one  D (shock)  boundary. 
Assuming  21  ~ 0 (~^-)  large  distance,  the  existence  of  the  last 
(double)  integral  requires 


<T  > 1 


(M.8) 


Is  is  assumed  tha t % lwr  vanishes  at  the  infinity;  this  is  satisfied 
by  (11.8). 


11.2  Remarks 


We  note  in  passing  that  for  the  airfoil  problem  governed  by  the  2-D 
transonic  sma  I I- d i s t urbance  equation,  we  have  £ = (i+<)  / Z K In 

this  case,  ^ = 0(  r ),  unless  the  circulation  is  zero;  thus;  2*  = 0 ( r'3 ) , 
/*-*•».  and  the  double  integral  exists.  At  large  r,  this  integral  can  be 
represented  by  a doublet 


~ Ifz**,*?. 


but  only  if  E i s integrable  ( C > 2)  , and  this  requirement  can  be  met 
only  by  a zero  lift.  In  general,  this  integral  is  proportional  to  xr‘7nr 
at  large  r (see  below). 

We  may  consider  two  functions  $ and  <PC  , satisfying  respectively, 
(with  z t £ ) 


v't-hZ', 


( I I .9.  11.10) 


which  may  not  fulfill  the  requirements  for  the  basic  relation  (11.7).  The 
latter  requirements  may  be  met,  however,  by  the  difference  if  S <f  - $c 
In  the  application  to  be  made  later  in  11.3  and  also  in  Appendix  III,  we 
will  choose  a function  £c  so  closely  approach  £ at  large  r that 


L-Lt-  O(f') 


<r  y 1 


(i 1 . 1 1) 


With  this,  we  construct  a > fulfilling  (11.2)  on  T and  (11.3)  on  D, 
also  approaching  at  large  r sufficiently  close  to  4>  , that 


v<p  = 7(*-<*t) 


as  r -+oo 


(11.12) 
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11.3  Captur i ,i,i  the  double  t^and  comparable  s i ngu  lar  itij-s 
i n the  tar  field  of  __  <pf 

The  for  field  behavior  of  <£#  (cf.  text  for  governing  equations) 
can  be  seen  from  the  relation  (11-7)  to  be  dominated  by  a vortex  of  a 
strength  JT0  = II$0I1TB  (assuming  a closed  2-D  body).  The  subsequent 
terms  of  orders  K'Ur  and  r'1  can  be  captured  with  the  help  of  $c 
and  (f  = <$oC  ■ with  2>*e  s(>  chosen  that  (T  > 2 . 

In  order  to  construct  the  anc*  Pc  ^or  t,,'s  case,  we  employ 

an  analytic  function  of  the  complex  variable  <*=xt  i'7 


Wcm<pc  + <PC  =/  Wc'<^ 
o 

iv' = 45  - c/p/t*  - ? - J 


(11.13) 


where  arj  (CJ-a)  and  arg.(^-  fc>)  are  restricted  to  the  range  (0,  Ztl^ 
caps  "a”  on  jro  . f , 0.  and  D have  been  omitted  for  convenience. 
This  function  Wc  has  the  desirable  properties:  (i)  -►  0.  as 

\$\-+0O  (ii)  W‘  is  continuous  across  the  boundary  T (T.V.  Sheet), 

(iii)  tlu-  circulation  around  the  slot  W (wing)  is  , and  the  total 

source  jA&rjU  d*  is  zero,  hence,  4>c  “►  l)  • as  I _*‘°°  ’ antJ 

(iv)  the**Kut  ta  condition  ( lfVc'  i co  ) is  satisfied  not  only  at  the  trailing 
edge  (x  = b)  , but  also  at  the  leading  edge  (V  = a).  We  then  have 

2.-  £»«  • •C|?r')-  o 

Now  construct  ^>aC  (omitting  "a”  from  <pol  ) from  f 


( I I . I A) 


which  admits  a particular  integral 


<£=  2 K'fo'-  f */»/] 

= _ liL  J i f (*,-*£)  f 

4K„  t 4 ; 


(II .15) 


*^Note  that^sy-i'^  and  that  is  the  partial  derivative^ 


withj  fixed, 
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The  4*oc  sought  is  simply 


rb 

+ J + fa\([+cB*]l"r'<lx,  (I  I -16a) 


which  is  analytic  at,  and  across,  the  entire  x-axis,  i.e.. 


I ^oc  H = H V hc/77  U = 0 


at  y = 0 , dll  x 


( II  . 1 6b) 


We  note  that  4G  is  continuous  at  Tf  = 0 and  that  source  distribution  in 
(11.16)  is  introduced  to  render  'Ofoe/?7}  continuous  by  cancelling  the 

jump  along  the  x-axis.  Thus,  we  can  apply  (II. 7)  to  <f>  ~d> 

“ ' ’0  *oC 

to  obtain  the  far-field  description  for  <pg  . Interestingly,  the  last 
integral  in  (I  I . 1 6 a ) does  not  give  rise  to  a source  in  the  far  field. 


This  is  because 


MM,)  ♦»«}, = o (11.17) 


on  account  of 


Kfe*  , j UclK(+c\  dx  = o 


i?0v'e  ?+  (V.  - we)  K ' ^ - (& )* • ^ 12 


+ z-hl £i)  (11.18) 

<p  /, 


the  final  result  confirming  (2.28)  in  the  text  is 


r f 


ft  * 5 Mf) * vt]  - N*l  ] * 


+ Q— . ♦ D/  — 

*/?|l  /$/* 


(II. 19a) 
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where  the  doublet  strengths,  however,  are  explicitly  related  to  the  solution 
interior  of  $ 

b 

Zir  D#l  = -jdhllxJx,  (M.I9b) 

o 

b 

2 n d;  - / K,  Vid  Z Bw  4*  - J d4<t>X-  4>' m \*x  + 

a a 

+2K”  (ll.19c) 

* 

,,  .i  , ,-*6 

We  point  out  that,  in  approaching  the  leading  edge,  (Pgx)  n-Const.  I X -A  | 

and  (Pen  ) \x-a\,  therefore,  the  last  (double)  integral  of  (Il.l9c) 

exist.  Observe,  also,  that  <^>c  and  its  derivatives  are  continuous  across 
D (the  shock),  and  (11.19)  is  unaffected  by  the  presence  of  the  (shock) 
d i scon  t i nu i t i es . 
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APPENDIX  III 


DETERMINATION  OF  THE  SOURCE  STRENGTH 

I I I . 1 Green ' s Theorem  Applied  to 
The  solution  $ = $o  * « satisfying  PDE  (2.6)  can  be  cast 

^ A . A 

forma  I I y in  a form  same  as  ( I I . 1 ) . In  terms  of  the  variables  £ - * *t 
and  <•  = k - i , and  droppi  ng  "a  ” f rom  , * , / , § and  % » (a  I so 

A 

writing  ? for  ).  The  PDE  in  question  is 

with 


<?,  AND  <3, 

$ = $o  ♦ 


The  third  term  of  arises  from  the  transformat  ion  of  (x,z,/)  to 

, y ) , with  K„  being  a function  of  ^ in  rj  -fKm  Z . In  the  following, 
we  will  apply  the  result  of  the  Green  Theorem,  (11.7).  to  the  function 

with  (poC  given  as  in  (11.16)  and  <P,c  so  constructed  to  fulfill  the 
requirement  for  applicability  of  (11.7).  Since  the  solutions  <po  and  £ 
have  only  the  weak  logarithmic  dependence  on  £ , the  result  obtained  below 
for  ( + € 4>t  ) can  be  equivalently  obtained  from  a separate  application 

of  (11.7)  to  ( <^>  - ) and  6 l $ ~ $(C  ) . The  reason  for  the  consideration 

of  their  difference,  y?  , is  for  the  economy  of  the  theoretical  effort,  since 
(11.7)  has  already  been  established  and  shown  to  be  unaffected  by  the 
presence  of  i mbedded  shocks  (for  Kn  ">  0 ) • The  app I i ca t i on  of  (11.7)  to 
e ( t - <P,c  ) for  the  case  in  which  the  shock  satisfies  the  transferred 
boundary  conditions  (2.26)  - (2.27),  does  not  appear  to  be  straight  forward. 
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Our  objective  is  to  capture  terms  in  the  far-field  description  for 
up  to?and  i nc I ud i ng f the  source  term.  For  this  purpose,  we  shall  find 
<p  among  the  solutions  to  the  Poisson  equation 


*4 


7C 


(I  I 1.3a) 


*ith  57  and  <t>  behaves  at  large  r as 
ic  lc 


zrz„  = o(i?n, 

be  = If  I ) . 


<r  > I , 


(II  I .3b) 

(I  I I -3c) 


For  convenience,^  and  its  ^-derivative  will  also  be  made  continuous 
at  the  x-axi s,  i .e. , 


I tie  B=  U = o , 


at  y)  =o  . 


(Ill -3d) 


The  key  equation  to  be  arrived  at  from  (11-7)  for  the  present  purpose  will 
be  therefore, 

,b 


(<Pe- <P0J*e(<Pr<t>.c')  = tr  - 

'a 


(III.1*) 


For  this  equation  to  be  valid  in  the  case  involving  a shock,  the  distribution 
of  Zj0  must  be  computed  with  the  discontinuities  of  Pox  and  <PqY^  evaluated 


(not  at  X — } but)  at  the  shock  boundary 


X=  e X^) 


(III. 5) 


Therefore,  the  relation  among  terms  with  subscripts  "0"  in  (ll 1.4)  is 
not  identical  to  the  corresponding  relation  for  the  leading  approximation 
obtained  previously  in  § 11.3,  Appendix  I I , as  long  as  g-*0.  In  the 
far  field  ( | -*ao  ),  however,  since  (cf.  § 11.3). 


*■  > 2 
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and,  whether  an  imbedded  shock  exists  or  not, 


Ifl 


Kl*  ■ 


(I  I I .6) 


Therefore,  the  presence  of  an  imbedded  shock  can  make  an  additional  contri 
bution,  if  any,  to  the  relation  in  the  far  field 


$-$,«  ~ c*" '($)*  2 *9"  - 

-rrffp  (ZrZ, . ' 0(1? n 


(III. 7) 


A 


only  as  a doub let , hence,  not  of  great  consequence.  Since  the  double 
integral  in  (III. 7)  vanishes  in  the  limit  |£  | — «.  00  when  condition 
(lll.3b)  is  met,  the  source  term  in  the  far-field  description  of  c 
must  be  the  same  as  in  that  of  <p  . The  remaining  task  is,  therefore, 
to  construct  Y and  <#>,  . 

Having  in  mind  the  requirement  (I  I 1.3b),  and  the  form  of  £ , we 

take  £ to  be 


> ©'  9 


OA 


'ith  and  chosen  to  yield,  at  large  | 

= O[(l"l?0ri?r*]i  y 


(III. 8) 


(I  I I .9) 


<P,-<P,cc=  OCO, 

where  p ^ od  . This  choice  of  Y will  make 

<-*i  c 

Z-  2„  = 


/ “It 

thus , fulfilling  ( I I I . 3b) 


(I  I I . 10) 


It  suffices  to  take 


or 


(111.11) 


where  <£c  and  <pgc  are  the  same  given  earlier  in  £ 11.3-  The  <pd  is  a 
Laplace  solution  (to  the  thin-airfoil  problem)  introduced  to  account  for 
the  doublet  behaviour  in  the  far-field  and  can  be  represented  by  the  vortex 
and  source  distribution  <pj  = R.P.  ( to^) 

b 

J -|TrJrf5  (HI.I2) 

a 

such  that 

f II  fax]}  d*  -H  = Oj  = 0 , 

Ja  » 

b 

j ([  fay  A xJx=  -TT  Dar 

a 


l in;. 


(II  I .13a, b, 


Thus,  at  large  /?/•  <PoA  * /?r2j  and  the  first  of  (ill -9)  is 

fulfilled.  Notice  that  the  "Kutta  condition"  at  the  leading  and  trailing 
edges  are  unnecessary  for  <pd  , as  long  as  (pj  x and  tj  are  integrable 
i n X.  on  the  wing. 

For  the  <f>  , i n Y , it  suffices  to  take  a particular  solution  to 

T / CC  t £ 

in  ( I I I . 8) , 


ncc  ^~>,c 

(III. 3a)  with  the  omission  of  the  nonlinear  (last)  term  of  £ 

- 1 


ic 


because  the  nonlinear  contribution  is  of  an  order  l£  I 
|£|  This  leads  to  (with  /<0=  K„  - (9  ) 


higher  at  large 


rf>  = 

T,cc  2K„ 


(in.14) 


With  <f>  and  (b  . so  chosen,  we  see  that  T,  can  fulfill  the 

T|  CC  '04  •—>  |C 

requirement  (111.10).  We  may  now  proceed  to  seek  a solution  to  (ill.  3a) 
for  the  <f>|C  . Thi  s w!  1 1 be  carr  i ed  out  i n I I I . 2 bel ow  for  stra i ght  obi  i que 
wing  only,  for  which  (£)'=  q 


c ,d) 
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I I I • 2 The*  far-field  description  and  determination  of  Q for 
straight  ob  I i gue  w i nqs 

With  the  foregoing  choices  of  <pgA  , 4>(  cc  . ',"'1  ^>c  2 can  l>e 


exprt 

have 


terms  of 

<t>c  ' 

2©  a , 

H-l 

irmfr**+ic; 

= Z®1 

±l(l 

K„  l 

Zdyl 

-WK”  - ( v-VcO*  -(*%)(#' 4 
+ ^ 0 M L I T * VC  - { */r  K - £*,"/’] 


4 *2*  ( W'J  w<"~  K'Wry-  Wc  Vfp 

(111.15) 

where  t lie  first  and  second  groups  of  terms  under  curly  brackets  arises  from 

(2©/K„)  3 fot/?  y ‘"ul  [ CHfO/Kr,  J (‘I’c  )x  ’ ( fi  f c ) x ’ respectively.  A 
particular  solution  <pf^  to  the  Poisson  equation  (I  I 1.3a)  can  be  obtained 
by  directly  integrating  with  respect  to  £ and  £ . adding  suitable  function 
of  ^ and  £ to  facilitate  fulfillment  of  conditions  (I  I 1.3c)  and  (I  I 1.3d) 
f°r  <pi  c.  The  results  can  be  written  as 

- *&§  I * (¥f,y(<^')  * !?>?)%( j*^)  - 

00  4) 


4 w/£Cfn/c')r-  u/c  - 
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(111.16) 


This  function  has  the  fol lowing  properties: 

(i)  At  large  |?J,  fulfills  not  only  (ill. 3c),  but  agrees  com- 

pletely with  4>  of  (2.29)  after  setting  0' = 0 therein,  except  for  the 
absence  of  the  (additive)  harmonic  functions  R.P  { _ <£Cr)  » £r*  ^ — ( Q + 

♦ --•*>}  ' 1 

(ii)  The  far-field  of  <J>jC  contains  no^  source  and  vortex  terms,  to  be 

sure . 

(P) 

(‘i’)  P/c  ,s  continuous  (no  jump)  across  the  x-axis  (i.e.,  at  7 = 0 
for  all  x)  , hence,  fulfilling  requirement  H <f>/c  ]]  = 0. 

We  observe  in  this  connection  that  ( WrWc)  Ind  ( Vj'c-  ^ ) are  synnetric 
in  y , so  are 


f M/eVe,  " Wc  ( §>iVc')y  - Wc(f  W'c). 

and 

*7  ("rV'Xtiq'-fWc)  i 


and  that 


can  be 


00 

represented  in  terms  of  ^ C * ) sff  JJ  = 4 & /f*-« Xx-»\\b  * 

‘W.  1 ^r~Wc— - 


l n gy  1 * jj.  (lvc-U/c)  [jjy (1<  (f  ip)) ! I f-tjdi 


which  again  has  the  symmetry  in  We>thus  may  complete  the  determination 

of  <J>{  by  adding  a source  d i s t r i but  i on  so  that 

= TTT  f Pv+.f/triE  1*  (111.17) 

*00 

may  fulfill  the  requirement  on  the  7>tt<  across  the  x-axis,  cf.  (III. 3d). 

In  approaching  the  x-axis,  ? can  be  evaluated  from 
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9 rx(rl\  - V3  \ ® i /./  ,1  x 

r,  (4' ) - 1 (Vf^)  “f\c  ~ 7.  »y  ( V 


V 


.<>  , 

~n(^Xxip  j 4cCs-,/)  S l»  /*-*)  <*$  \ (,|,.l8) 

/l 


This  gives 


fpl 


^,1=  ?rf 

,b 

-%  Wct'Iltjj  1C(S  „)S  l»lx-s/Js 


(111.19) 


We  observe  that,  over  the  wake  portion, 

.(P\ 


ff Hc/?7 H=  o t 


(III .20) 


since  ^ and  cf>cx  is  continuous  across  the  x-axis  downstream  of  the  T.E. 

Since  c^c  does  not  give  rise  to  a source  in  the  far-field,  we  have, 

as  l?|-*  00  . a source  strength  for  c p as  (recalling  "A"  has  been 
/£  a . ,c 

omitted  from  <p  and Q ) 


/ /») 

Q,*3,£=-*irj  /*n  JJ  d* 


= 2tt  © a m /» « ±rl  «. 

/<.  3/  *+  8*rKlUf7  »-a 
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4 


* 77*  ® l f “Ml  W * 

ZK”  J<x 


8K„zrr  <*/  4^  A 7 


Jj± 

anK1, 


b t> 

\ ©/ (Lkxxfl^  j 1e(S>,)S  l„lK-sldl 

M /)  ' 4 


(111.21) 


The  above  result  for  straight  oblique  wings(with  ©•=o)  does  not 
require  the  assumption  of  a local  similarity  in  the  3-D  structure.  In 
cases  in  which  the  3_D  similitude  applies,  the  far-field  source  strength 
CJ(  for  the  similarity  solution^,  cf.  ( 4 . 1 4b)  , can  be  determined  f rom 
( I I I . 2)  through  (4 . 1 6b) 


'N^  (3 

0 — - — - -t 

' ©CC'  4 7TK 


- IT,  Uc 

Z • 


noting  in  this  case  that  $'/£) I**/* ? and  using  (4.1 5) 

to  relate  fT  , D , and  £)1  to  [f'  , Qr  , and  D * ■ The  final  result  for 

11 0 0 o ’to  o 0 


the  source  strength  in  the  far-field  £>  is  (with^-a^  2) 


Q.  = 


- ^ [3, 


K. 


sttk: 


(111.22) 
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FIGURE  CAPTIONS 


Figure  1 Illustration  of  the  reference  chord,  the  span  and  the  aspect 
ratio  used  for  a straight  oblique  wing(a),  a more  general  curved  ob- 
lique wing(b),  and  symmetric  swept  wings(c) ; the  shaded  areas  in  (c) 
indicates  regions  of  breakdown.  Breakdown  (nonuniformity)  of  the  theory 
may  also  occur  in  the  tip  region  (cf.  text). 


Figure  2 Illustration  of  the  Cartesian  and  the  orthogonal  curvilinear 
coordinate  systems  in  the  wing  plane  ( ZzZ  z 0 ) ■ 

Figure  3 Examples  of  *o*  ’ <Pl  K , and  tzi  on  the  upper  and  lower  surfaces 
of  a straight  oblique  wing  at  3 • 6 . The  airfoil  section  is  NASA 
3612-02,  40  rescaled  to  an  arbitrary  thickness .(  Refer  to  text  for  results 
shown  in  dots);  wing  has  no  twist  and  bend;  local  incidence  is  zero. 


Figure  4(a)  Results  for  a straight  oblique  wing  at  «„=  3-45  at  two  inci- 
dences : £ ot  = 0 (in  solid  curve)  and  a OUs.  0 . 532  T ( i n dash  curve). 

Distribution  of  *.*• 

Figure  4(b)  (Continued)  Distribution  of  t\ x . 

Figure  4(c)  (Continued)  Distribution  of  • 

Figure  5(a)  Comparison  with  a major  to  minor  axes  ratio  of  20  at  22.5° 

sweep  and  a free-stream  Mach  number  0.8242.  The  present  result  is  shown 
in  solid  curve.  This  graph  g i ves  C_  d i s t rubut ion  at  the  span  station 

7 =.-0.69.  P 

Figure  5(b)  (Cont  inued)  Span  station  y =»  0. 


Figure  5(c)  (Cont i nued) Span  station  y"=.69- 


Figure  6(a)  Comparison  for  a 
minor  axes  ratio  of  14  at 
surface  Cp  at  span  station 

Figure  6(b)  (Continued)  Span 

Figure  6(c)  (Continued)  Span 

Figure  6(d)  (Continued)  Span 

Figure  6(e)  (Continued)  Span 


12%  thick  elliptic  oblique  wing  with  major  to- 
30°  yaw  and  = 0.7549-  This  graph  gives 

7 = - 0.89. 

stat ion  y s -O.69. 
stat  ion  y =0. 
station  y=0.69. 
station  y = 0 . 89 


Figure  7(a)  Comparison  for  a 12%  thick,  elliptic  oblique  wing  with  an  axes 
ratio  1 4 at  30°  yaw  and  = 0.7677  This  graph  gives  surface  C«  at  span 
station  y r-0.79- 
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F i gure 

7(b) 

(Cont i nued) 

Span 

station 

Figure 

7(c) 

(Cont i nued) 

Span 

station 

Figure 

7(d) 

(Cont i nued) 

Span 

station 

Figure 

7(e) 

(Conti nued ) 

Span 

station 

F i gure 

7(f) 

(Cont i nued) 

Span 

s tat  ion 

Figure 

7(g) 

(Conti nued) 

Span 

station 

Figure  8 Surface  distributions  of 

an  airfoil  section  NACA  6^)006  for 

K„=  2.50.  The  surface  value  of 
distribution  shown  (cf.  text). 


7=  -0.59- 
1 - -0.20. 
y=  0. 

7 = 0.20. 
y ao.50. 
y = o.79. 

ftjr  on  a nonlifting  oblique  wing  with 
a component  similarity  parameter  of 
can  be  recovered  from  the  Cp'^ 


Figure  9 Compar i son  of  surface  on  a nonlifting  elliptic  oblique  wing 
obtained  from  loca I -s imi lar i ty  solutions  (in  discrete  symbols)  with 
corresponding  numerical  data  via  the  unsteady  analogy  (in  solid  curves). 
The  airfoil  section  and  K*value  are  the  same  as  in  Figure  8,  and  £0  is 
taken  to  be  0.10,  using  ^>,5-  for  case  (A)  of  Figure  8.  The  Cp  is  computed 
for  a 6%  thickness  ratio  and  M„~  0.827-  The  latter  conditions  together 
wi  th  €0  = 0.10,  gives  VI  = 0.895  and  ^ = 22.6  for  a sweep  angle  of 

22. 5°.  * ' 

Figure  1 0 ( a) Compar i son  of  surface  Cp  on  a nonlifting  obique  wing  of  len- 
ticular planform  from  the  1 oca  I -s imi 1 a r i ty  solution's  (in  discrete 
symbols)  with  the  corresponding  numerical  data  via  the  unsteady  analogy 
(in  solid  curves).  The  airfoil  section  and  K*  value  are  the  same  as  in 
Figure  8,  and  €0  =0.10,  using  ft?  for  case(A)  of  Figure  8.  This  graph 
shows  results  for  three  span  stations  on  the  aft  panel,  y"  = 0,  0.5  and 
0.8. 


Figure  10(b)  (Continued)  Span  stations  on  the  fore  panel , y*  = -0 . 2 and  -0.5. 

Figure  1 1 (a) Compar i son  of  surfaceCp  on  a nonlifting  oblique  wing  of 
ha  1 f- 1 ent icul ar  planform,  with  the  straight  axis  at  the  leading  edge. 

The  ft?  for  case  (B)  is  used  for  computing  Cp  which  is  computed  other- 
wise in  the  same  manner  as  in  Figure  10.  This  graph  shows  results 
for  span  station^z  0,  0.50  and  0.80  on  the  aft  panel. 

Figure  1 1 (b)  (Cont  inued)  Span  s ta t ions  y’ -0 . 2 , and  -0.50  on  the  fore  panel. 

Figure  1 2 (a ) Compar i son  of  surfaceCp  on  a nonlifting  oblique  wing  of  half- 
lenticular  planform,  with  the  straight  axis  at  the  trailing  edge.  This 
pertains  to  case  (C)  in  Figure  10.  The  graph  shows  results  for  span 
stations  p's  0,  0.50,  and  0.80. 

Figure  1 2 (b) (Cont i nued)  Span  stations  y"  - -0.2  and  -0.5  on  the  fore  panel. 

Figure  11.1  Illustration  of  the  wing  boundary  (W) , the  trace  boundary  (T) , 
and  the  shock-discontinuity  boundary  (D)  . 
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and  mixed,  being  governed  by  equations  consistent  with  the  transonic  small- 
disturbance  approx imat ion . The  work  analyzes  3-D  corrections  to  this 
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‘nonlinear  problem  and  matching  its  solutions  to  that  of  an  outer  flow.  In 
the  (parameter)  domain  of  interest,  the  outer  solutions  correspond  to  a 
high  subsonic,  or  a linear  sonic,  outer  flow,  representable  by  a 
Prandt l-Glauert  solution  involvinn  a swept  (or  curved)  lifting  line  in  the 
leading  approximation. 

Among  the  3'D  effects  is  one  arising  from  the  compressibility  correc- 
tion to  the  velocity  divergence,  absent  in  classical  works.  More  important 
are  the  upwash  corrections  which  include  the  influence  of  both  the  near  and 
far  wakes,  as  well  as  the  local  curvature  of  the  center  line.  For  straight 
oblique  wings,  local  similarities  exist  in  the  3"0  flow  structure,  permitting 
the  reduced  equations  to  be  solved  once  for  all  stations.  An  analogy  also 
exists  between  the  oblique-wing  problem  and  that  of  a 2-0  transonic  flow 
which  is  weakly  time-dependent;  this  provides  an  alternative  method  of 
solving  numerically  the  inner  airfoil  problem.  Of  theoretical  interest  is 
the  uncovering  of  a compressibility  correction  to  the  outer  flow  in  the  form 
of  a line  source  and  a nonuniformity  of  the  perturbation  solutions  at  a 
shock  root  (and  its  treatment). 

A procedure  based  on  a line  relaxation  method  for  solving  numerically 
the  reduced  inner  problem  is  described;  solutions  with  high  subcritical,  as 
well  as  slightly  supercritical,  component  flows  are  demons t ra ted . Comparison 
with  corresponding  numerical  solutions  based  on  f u 1 I -potent i a I equations 
for  oblique  elliptic  wing  shows  encouraging  agreement. 
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